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EXECUTIVE SUMMARY 




I INTRODUCTION 


A 1979 report by the National Aeronautics and Space 
Administration's Atmospheric Lidar Working Group defined seven scien- 
tific objectives that could be usefully addressed by an orbiting lidar. 
The objectives are: 

(1) Defining the global flow of water vapor and pollutants. 

(2) Improving stratospheric and mesospheric chemistry and 
transport models. 

(3) Improving atmospheric radiation-balance models. 

(4) Augmenting the meteorological data base. 

(5) Understanding upper atmospheric waves. 

(6) Investigating thermospheric atomic species. 

(7) Investigating magneto spheric aspects of sun/weather 
relationships • 

Aerosol and cloud measurements play vital roles in Objectives (l)-(3) 
and (5), and an important, though secondary, role in Objective (4). In 
addition, stratospheric/ tropospheric temperature profile measurements 
are important to each of Objectives ( 1 )— (5). 

The Working Group report* defines the accuracy and resolution that 
measurements must have for each objective. In addition, the report 
briefly indicates the accuracy and resolution thought to be achievable 
in various orbiting lidar experiments. However, the Working Group's 
activities did not include detailed and complete calculations of aerosol 
and cloud retrievals for the wide variety of atmospheric and measurement 
conditions that are likely to be encountered in practice. Thus, the 
question was not fully explored as to how well shuttle lidar aerosol and 
cloud measurements could contribute to each objective. Also not 
explored was the possible use of lidar backscatter profiles at the tri- 
pled Nd:YAG laser wavelength (0.355 ym) to aid aerosol measurements and 
possibly to provide temperature profile measurements. 
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During the Working Group meetings. General Electric Space Division, 

o 

under contract to NASA, began a facility definition study 4, to better 
define the lidar system hardware that would be compatible with both the 
space shuttle accommodations and the Working Group's scientific objec- 
tives. That study established hardware parameters that could be used in 
scientific simulation studies with the confidence that they could actu- 
ally be implemented in a future shuttle-borne system. 

The Atmospheric Lidar Working Group report and G.E. Facility Defin- 
ition report thus laid the groundwork for the present study. The objec- 
tive of this study is to assess quantitatively, for a variety of realis- 
tic scenarios, the accuracy and resolution with which aerosol, cloud, 
density, and temperature profiles could be retrieved from measurements 
made by an orbiting Nd:YAG lidar having parameters consistent with the 
G.E. study. This process indirectly identifies the lidar parameters 
required to obtain scientifically useful data. The aerosol and cloud 
assessments were to be made using an error-analysis-and-simulation code 
available at SRI; the code was to be improved as part of the study. The 
density and temperature assessments were to be made using a new analysis 
technique based on elastic backscatter at wavelength 0.355 ym (with 
aerosol and cloud corrections from longer wavelengths). This technique 
yields density profiles that are then integrated vertically to yield 
pressure and temperature profiles. Subsidiary goals were to: 

• Define a realistic set of aerosol, cloud, density, and ozone 
models for input to the assessment. 

• Define the correlative sensors that could best aid the lidar in 
achieving the objectives of the Working Group report. 

• Point out specific hardware implications of the necessary lidar 
parameters . 

• Expose the study methods and results to the scrutiny of the 
scientific community, for example by publishing them in a 
refereed journal. 

A specific question addressed by the study was whether the scientific 
benefit expected from the measurements at 0.355 ym would warrant the 
additional complexity of transmitting, detecting, and processing signals 
at that wavelength. 
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Part 2 of this report gives detailed results of the aerosol and 
cloud analyses, while Part 3 gives the results of the density and tem- 
perature (plus improved aerosol and cloud) analyses. Those parts are 
formatted for rapid conversion to journal manuscripts, in support of the 
last objective mentioned above. The remainder of Part 1 summarizes the 
results and makes the recommendations mentioned above. 

This report has benefitted considerably from discussions with 
E.V, Browell and S.T. Shipley, who also performed a critical review of a 
previous draft version. 
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IX 


RESULTS 


A* Atmospheric Models 

Atmospheric models were developed for low, middle, and high lati- 
tude bands. The models include profiles of total gas density, ozone, 
cirrus and noctilucent clouds, as well as mesospheric, stratospheric, 
and tropospheric (Saharan and marine) aerosols. Profiles of back- 
scattering and extinction coefficient at wavelengths 0.266, 0.355, 

0.532, 0.694, and 1.064 ym were computed for each model. The models are 
based on previous optical and physical measurements, and have been com- 
pared to previous ruby lidar measurements, so as realistically to 
represent the conditions encountered in practice. The models and their 
derivation are described in Appendix A of Part 2. 


B. Lidar Parameters 

A set of lidar parameters was developed after careful review of the 

2 

G.E. Facility Definition report, the scientific measurement goals, and 
our own review of hardware capabilities. The detailed parameters are 
listed in Table 1 of Part 2. Important points to note here are the 
following : 

• Three wavelengths (1.064, 0.532, and 0.355 ym) are transmitted 
simultaneously, but they are obtained from a single Nd:YAG laser 
by harmonic generation, with no special line-narrowing, tuning, 
or wavelength-stability requirements. 

• Three detectors and two dichroic beamsplitters are used to mea- 
sure all three wavelengths simultaneously, but detector 
linewidth and stability can be provided by interference filters, 
with no need for interferometers or other more sophisticated 
spectral discriminators. 
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• Transmitted energies and beam divergences satisfy the prelim- 
inary eye safety guidelines* developed by G.E. 

• The weight, power, space, and other requirements implied by the 
lidar parameters are consistent with the limits for shuttle sci- 
ence payloads. 

In sum, the lidar system assumed here may be considered as a baseline 

system, suitable for an early shuttle mission. The transmitter is less 

sophisticated than the Nd :YAG-pumped dye laser system considered in the 

i o 

Working Group and Facility Definition reports * as the first module in 

an evolutionary system. Nevertheless, the engineering aspects of 

transmitting the fundamental, doubled, and tripled wavelengths simul- 

o 

taneously from shuttle have not been worked out. Also, the three- 

wavelength detection system is more complex than the two-wavelength 

2 

systems emphasized by the Facility Definition report. Finally, we have 
followed the Facility Definition report in assuming different receiver 
fields of view and filter bandwidths for viewing daylite and dark por- 
tions of the earth. (See Table 1, Section IV-B of Part 2). This pro- 
cedure is necessary to maximize scientific gain consistent with eye 
2 

safety; however, it implies the use of a system that automatically 
changes telescope beam stops and interference filters upon crossing from 
dark to daylite areas* 


C. Other Inputs 

Background radiances at the three lidar wavelengths were developed 

to represent sunlit clouds, sunlit ocean, and moonlit clouds, each as a 

function of source (sun or moon) zenith angle. (See Table 2, Section 

IV-C of Part 2.) Probable errors in atmospheric transmission and in 

calibrating the received lidar signal in terms of atmospheric back- 

3 

scatter were adapted from a previous study, taking into account the 
larger range of heights and aerosol and cloud types considered in the 
current study. In particular, we adopted a lcr uncertainty for one-way 
transmission in an aerosol or cloud layer of one-half the particulate 


*See also Section IV-B of Part 2. 
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optical thickness x p, because effective backscatter-to-extinction ratios 
in most appreciably turbid media (i.e., those with >_ 0.3) cannot be 
estimated to an accuracy of better than ±50 percent, even with the aid 
of multiwavelength lidar measurements. (See Appendix B of Part 2.) We 
also assumed that the root-mean-square (rms) errors in conventional 
molecular density profiles (used to analyze the lidar data) were 
appropriate to data-sparse regions of the globe where most shuttle lidar 
data will be acquired and where shuttle lidar is most likely to fill 
data gaps. Specifically, we assumed conventional density rms errors 
ranging from 2 to 3 percent below 30 km and 5 to 10 percent above 40 km. 


D. Aerosol and Cloud Measurements at 0.53 and 1.06 pm 

We first considered a lidar with only the fundamental and doubled 
wavelengths (1.06 and 0.53 ym) available, each being analyzed indepen- 
dently, using conventional density data, to retrieve aerosol and cloud 
backscatter profiles. Simulations were conducted for low, middle, and 
high latitude model atmospheres, assuming a variety of daytime and 
nighttime backgrounds. For each scenario, expected errors in retrieved 
particle backscatter were calculated from algebraic expressions and also 
checked by numerical simulations using random number generators. Both 
the algebraic expressions and the simulations include the four sources 
of error that we have encountered in actual lidar measurements and 
analyses — namely, errors in (1) measuring the backscattered signal, 

(2) accounting for transmission losses, (3) determining the total gas 
density profile at the lidar location, and (4) calibrating signals in 
terms of atmospheric backscatter by using an atmospheric layer where gas 
backscatter dominates particle backscatter. 

The results show that useful retrievals can be made for the follow- 
ing constituents and conditions: 

• By day : Vertical structure of tenuous clouds (subvisible and 

visible), Saharan aerosols, and boundary-layer aerosols (at both 
0.53 and 1.06 ym wavelengths), and strong volcanic stratospheric 
aerosols (at 0.53 ym) . Quantitative backscatter can be 
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determined for all of the preceding, provided particulate opti- 
cal depth does not exceed -0.3. 

• n ight : Vertical structure and quantitative backscatter for 

all of the above, plus upper tropospheric and stratospheric 
aerosols (at 1.06 pm) and mesospheric aerosols (at 0.53 pm) and 
noctilucent clouds (at 0.53 and 1.06 pm). (Again, quantitative 
backscatter retrievals suffer when particulate optical depth 
exceeds -0.3). 

These results were obtained with vertical resolutions of 0.1 to 0.5 km 
in the troposphere and 0.5 to 2.0 km in the stratosphere and above, 
except within the mesospheric aerosol and noctilucent cloud layers, 
where resolutions of 1.0 and 0.25 km, respectively, were used. Horizon- 
tal resolution was 100 km in the troposphere, 200 km in the stratosphere 
and lower mesosphere (below -60 km) , and 2000 km in the upper mesosphere 
(to -85 km). Vertical structure of the stronger layers (noctilucent and 
cirrus clouds, lower tropospheric aerosols) could be observed on much 
finer horizontal scales (sometimes on single shots — about 700 m apart); 
moreover, these fine-scale results could be converted to quantitative 
backscatter if a transmit energy monitor or stable transmit energy were 
pr ovid ed . 

In general, the principal virtue of the 1.06-pm wavelength is its 
reduced sensitivity to gas density errors, which are the dominant error 
source for upper tropospheric and nonvolcanic stratospheric aerosol 
retrievals. Also, 1.06-ym retrievals are somewhat less susceptible to 
transmission errors in the stronger layers. The principal virtue of the 
0.53-ym wavelength is the stronger signals available (because of atmos- 
pheric backscatter and quantum efficiency). These stronger signals are 
essential for observing the mesospheric aerosol layers at useful hor- 
izontal and vertical resolutions. 

At 0.53 and 1.06 ym, errors in the gas density profile are the dom- 
inant source of uncertainty in measuring nonvolcanic stratospheric and 
upper tropospheric aerosols. Although useful retrievals of these aero- 
sols can be made at 1.06 ym using conventional density information, 
retrievals at both wavelengths would be Improved by better density 
information. 
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E. Density , Temperature , Aerosol , and Cloud Measurements 

at 0. 355 , 0. 532 , and 1. 064 ym 

We next considered a lidar having three wavelengths — 0.355, 0.53, 
and 1.06 ym. An analysis technique (Figure 1 of Part 3) was developed 
that combines information from 0.355 and 1.06 ym to yield a relative 
density profile and an uncertainty profile that takes into account aero- 
sol and cloud contamination. The relative density profile is converted 
to an absolute density profile by normalization at an assumed isopycnic 
level or a level where conventional density data are available to good 
accuracy. The lidar-derived density profile is used to improve aerosol 
and cloud retrievals at 0.53 and 1.06 ym, and is also integrated verti- 
cally (from a reference height where temperature or pressure are 
guessed) to yield pressure and temperature profiles. The two keys to 
the success of this analysis procedure are: 

• Returns at 0.355 ym are much more sensitive to gas backscatter- 
ing than are those at 1.06 ym. (Hence, on a relative basis, 
0.355 ym returns are much less sensitive to aerosol contamina- 
tion, whereas 1.06 ym returns readily indicate regions where 
particulate contamination might be significant.) 

• Aerosol retrievals and the lidar-inferred temperature profiles 
are relatively insensitive to bias errors in the density pro- 
file, so that errors in absolute density normalization have very 
little effect on these data products. Also, this reduces the 
effect of some transmission errors. 

Error analysis equations were developed to quantify the above argu- 
ments and to provide error bars on all retrieved quantities. Also, to 
check the algebraic error expressions, the numerical simulation code was 
extended to inject appropriate errors into the multiwavelength retrieval 
procedure. In addition to the four error sources listed in Section II-B 
(signal, transmission, conventional density, and calibration), we also 
included errors in absolute density normalization, in the reference tem- 
perature or pressure, and in estimating short-wavelength particle back- 
scatter from a long-wavelength measurement (with errors). 

Simulations were run for the low, middle, and high-latitude model 
atmospheres. In general, cirrus clouds (subvisible as well as visible) 
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and lower tropospheric aerosols introduced density and temperature 
errors that were larger than those in conventional data; daytime signal 
measurement errors were similarly large in the stratosphere and above. 
Hence, useful lidar density and temperature measurements were restricted 
to nighttime conditions in the stratosphere and (cloud-free) upper tro- 
posphere. (Signal errors also prevented useful measurements in the 
mesosphere if horizontal resolution was 2000 km or less.) 

Simulations were run with vertical resolutions of 0.5 to 1.0 km in 
the upper troposphere and lower stratosphere, and 2.0 km in the upper 
stratosphere. A horizontal resolution of 200 km provided useful signal 
accuracies (0.5 to 2.0 percent) below 40 km; however, this was increased 
to 2000 km to extend useful signals to 55 km. Retrieved relative den- 
sity profiles had rms errors of 0.5 to 2 percent in the upper tropo- 
sphere and stratosphere. The limiting error at most heights was signal 
measurement error. However, nonvolcanic stratospheric aerosols intro- 
duced errors of 1 to 2 percent at the mixing ratio peak (typically _27 
km at low latitudes to ~16 km at high latitudes). For a given amount of 
particle backscatter, the density error tends to be larger at low lati- 
tudes, because the low-latitude peak occurs at greater heights, where 
the smaller absolute density yields larger relative errors. A midlati- 
tude moderate volcanic stratospheric aerosol model introduced relative 
density errors of about 3 percent at the peak (~18 km). 

Use of the lidar-derived relative density profiles in place of con- 
ventional profiles significantly improved the accuracy of aerosol 
retrievals in the nonvolcanic stratosphere and upper troposphere. This 
was especially so at 0.53 pm, where conventional density errors had led 
to very poor retrievals at these heights. The impact of lidar-derived 
density errors on volcanic stratospheric aerosol retrievals was minimal, 
both because of the insensitivity of such retrievals to density errors 
and because of the poorer accuracy of the volcanic lidar density 
profile • 
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For cloud-free, nonvolcanic conditions, lidar-derived temperature 
profiles had ms errors of 1.2 to 2.5 K in a layer bounded on the bottom 
by strongly scattering tropospheric aerosol layers and on the top by a 
height ~8 km below the reference height where temperature or pressure 
was guessed. For the model atmospheres used here, this lower bound was 
~5 km; the upper bound was ~32 km for a horizontal resolution of 200 km 
and ~47 km for a horizontal resolution of 2000 km (i.e., the reference 
heights were at ~40 and ~55 km, respectively). The midlatitude moderate 
volcanic stratospheric aerosol introduced temperature errors of ~3 K 
between 17 and 20 km. These simulations assumed an 8 K ms error in the 
temperature guessed at the reference height (40 to 55 km). 

The major advantage of these lidar-derived temperature profiles 
over those obtained by passive nadir-viewing spacebome sensors is their 
fine vertical resolution (-0.5 km in the upper troposphere and tropo- 
pause region, increasing to -2 km in the upper stratosphere). This 
resolution would permit defining the tropopause and temperature-wave 
structures to a degree of detail never before possible. 

In general, the simulations validated the algebraic expressions 
used to put error bars on lidar-derived density and temperature pro- 
files. This is important, because the errors reported above apply to 
specific atmospheric situations and will, for example, increase consid- 
erably whenever strong particulate (cloud or aerosol) layers are encoun- 
tered. In this respect the outstanding feature of the multiwavelength 
analysis technique is that it simultaneously retrieves density and par- 
ticle profiles, and includes particle contamination effects in the den- 
sity and temperature error bars. Thus, the user is immediately aware of 
significant cloud or aerosol contamination. 
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Ill RECOMMENDATIONS 


A. Hardware Development 


A major conclusion of this study is that considerable scientific 
benefits can be obtained by adding measurements at 0.355 ym. Therefore, 
we recommend that: 

• A ground-based or airborne three-wavelength NdsYAG lidar system 
be built to make aerosol, cloud, density, and temperature mea- 
surements, both to test the conclusions reached in this study 
and to aid in refining the analysis techniques by using actual 
data. 

• The three-wavelength Nd:YAG transmitter/detector configuration 
be included in any future engineering studies of the shuttle 
lidar facility. 


This study has also shown that, in the lowest few kilometers of the 
atmosphere, aerosol structure can be usefully retrieved with a vertical 
resolution of 0.1 km or less, sometimes from single-shot returns. Else- 
where in the troposphere, even subvisible cloud heights can be retrieved 
to 0.25 km or less. Vertical resolution of 0.5 to 1 km is useful in the 
lower stratosphere, and 2 km above, with 1 and 0.25 km vertical resolu- 
tion in the mesospheric aerosol and noctilucent cloud layers, respec- 
tively. This implies that: 


• The shuttle lidar data recording and/or telemetry system be 

designed to preserve three simultaneous profiles for each shot, 
with vertical resolution similar to the following: 


Altitude 

(km) 

0-5 

5-18 

18-30 

30-40 

40-75 

75-85 

85-90 


Vertical 

Resolution 

(km) 

0.1 

0.25 

1 

2 

1 

0.25 

2 
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Most signal profiles will be partly in the pulse- counting regime, 
partly in the current regime, and partly in the intermediate regime (-1 
to 100 pulses p, s 1) where both of these conventional recording tech- 
niques encounter practical difficulties (e.g., pulse overlap and nonuni- 
form charge per pulse). The data analysis techniques require that each 
part of any profile be relatable to other parts with an accuracy of 
-0.5% or better. Therefore we recommend that: 

• A reliable, repeatable, and precise method^ of logging and 

splicing together all parts of the signal profile from 0 to 90 
km be developed and tested. The method should permit determina- 
tion of the ratio of signal return from any height to that from 
any other height to an accuracy of 0.5% or better, as well as 
accurate logging of the strong surface return. (The above accu- 
racy excludes the jlT statistical error inherently determined by 
the pulse-generation rate in the pulse-counting mode; this error 
can be reduced by summing many shots or expanding the size of 
pulse-counting range bins.) 


B. Further Studies 


This study has raised several important questions that could be 
answered by further study. Therefore we recommend that: 

• The simulation program be used to evaluate the expected perfor- 
mance of: 

- A ground-based or airborne three-wavelength lidar design, as 
recommended in Section A. 

- A two-wavelength lidar based on the alexandrite solid-state 
laser (A =0.7 to 0.8 pm) and second harmonic generation. 

Such a system would benefit from increased quantum efficiency 
and atmospheric backseat ter at -0.75 ym as compared to 1.06 
pm, while retaining many of the benefits of the more complex 
three-wavelength system simulated in this study. The tunabil- 
ity of the alexandrite laser also favors other measurements, 
such as water vapor by differential absorption. 

- A three-wavelength spacebome lidar with increased output at 
the (relatively eye-safe) 0.355 ym wavelength, to extend the 
upper height limit and improve the horizontal resolution of 
density and temperature measurements. 

• A statistical study of existing radiosonde data be conducted to 
determine the rms error with which density can be estimated at 
the 8-km "isopycnic" layer, as a function of latitude and 
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season. Also, that existing and new data (from dust sondes, 
lidar, and satellites) be surveyed to estimate how often signi- 
ficant aerosol and cloud contamination occurs at this level. 

• A study be conducted to explore the use of three— wavelength 
lidar data (with conventional density data) in the lower tropo- 
sphere to estimate aerosol size distributions and backscatter- 
to-mass conversion factors. 

• A study be conducted to determine whether the wavelength- 
dependence ratio ¥ [= B (0.355 pm) /B (1.064 pm)] can be 
estimated to better than the ±100% accuracy assumed in this 
study by using 

- Height information together with previous data on the expected 
height-dependence of ¥ • 

- Information on Bp(0.532 pm). 

• The existing simulation code be improved to 

- Make input parameters easier for the user to locate and 
understand . 

- Handle random errors that produce a negative inferred signal. 

- Include effects of NO 2 absorption. 


C. Correlative Sensors 

We recommend that the following correlative sensors be operated in 
conjunction with the shuttle lidar, both to augment the scientific data 
set and to validate the lidar data: 

• Ojl fc h e Shuttle : 

— A radiometer of the Earth Radiation Budget type, to measure 
radiative effects of aerosol and cloud layers detected by the 
lidar. 

- A cloud physics radiometer,^ to compare passively inferred 
cloud properties (height, phase, temperature, optical thick- 
ness) with those inferred by lidar. 

- A temperature profiling radiometer, to permit comparisons of 
lidar and passively inferred temperature profiles, and studies 
of aerosol and cloud effects on each. 

- One or more surface property sensors, such as a Coastal Zone 
Color Scanner or a Landsat Mapper, to study effects of aero- 
sols and thin clouds on measured radiances and inferences. 

• Elsewhere : 

- Radiosondes, to measure temperature and density profiles for 
comparison to model, interpolated, and lidar-derived profiles. 
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Dustsondes, to measure stratospheric and upper tropospheric 
aerosol profiles. 

Ground-based and airborne lidars, to measure clouds, density, 
and temperature. 

In situ aerosol samplers (impactors, filters), to determine 
aerosol optical microproperties (size distribution, composi- 
tion, shape). 

Satellite limb-scanning aerosol radiometers (such as SAGE II), 
to measure aerosol extinction profiles. 
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Part 2 

ORBITING LIDAR SIMULATIONS: 
AEROSOL AND CLOUD MEASUREMENTS 
BY AN INDEPENDENT -WAVELENGTH TECHNIQUE 




I INTRODUCTION 


Previous studies have considered the use of an orbiting lidar to 
make aerosol and cloud measurements. Several of these studies ° have 
estimated the measurement accuracy of lidar return signals for various 
orbital altitudes, lidar parameters, and atmospheric conditions and 
illuminations. However, to date, no quantitative studies have shown the 
accuracy with which cloud and aerosol backscattering profiles, and other 
parameters such as extinction and mass, could be retrieved from the 
lidar signal profiles for a variety of realistic scenarios. 

The difference between signal-measurement accuracy and particle- 
retrieval accuracy is often crucial. For example, typical lidar parame- 
ters and atmospheric conditions often yield lidar signal profiles that 
are more accurate at midvisible wavelengths than in the near-infrared; 
nevertheless, errors in gas backscattering and atmospheric transmission 
often yield particle backscatter information that is far less accurate 
in the midvisible than the near-IR, in spite of the more accurate mid- 
visible signals. Other degradations are also possible. Because studies 
of aerosol and cloud formation, transport, and radiative and climatic 
effects will use the retrieved particle profiles, rather than lidar sig- 
nal profiles, it is important to estimate the quality of these retrieved 
profiles in judging the utility of shuttle lidar measurements. 

This paper presents the results of simulated lidar measurements and 
retrievals for several realistic combinations of lidar parameters, aero- 
sol and cloud profiles, background lighting conditions, and retrieval 
error sources. These results not only quantify the expected errors in 
retrieved data, but also highlight areas that require early hardware 
development or further study. 
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We refer to the simulations presented here as "independent- 
wavelength" simulations. This terminology describes retrievals in which 
the lidar data are analyzed independently at each wavelength, even 
though other lidar wavelengths might he producing data simultaneously. 

An analysis technique that combines signal profiles from two or more 
wavelengths to retrieve gas density, temperature, aerosol, and cloud 
information is described and simulated in a companion paper. ^ In both 
this and the companion paper we have taken care to make the simulated 
analysis techniques as similar as possible to techniques used with 
actual lidar data in the past, and to include error sources that real- 
istically mimic those encountered in actual measurements and retrievals. 
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II SUMMARY OF SINGLE -WAVELENGTH 
AEROSOL LIDAR ERROR ANALYSIS METHOD 


The error analysis method used here and described in a previous 

paper is based on the scattering-ratio/normalization method of lidar 

calibration and data analysis. This method has been used, in varying 

8—20 

forms, in many lidar studies. Here, we summarize the analysis method 
and the resulting algebraic error expressions. 

A fundamental quantity that occurs in this analysis method is the 
scattering ratio, defined by 

B p ( X , z) + B g (A,z) B p (A ,z ) 

R(A,z) = 1 + , 

Bg(A ,z) Bg(.\,z) (1) 


where B p and Bg are the particulate and gaseous backscattering coeffi- 
cients, respectively, at wavelength X and height z. As shown in detail 
elsewhere, the scattering ratio is related to backscattered lidar sig- 
nal S(A,z) as 

2 

(z - z L ) S(A,z) 

R(X,z) = , 

K(A)Q (A ,z,z l )D(z) (2) 


where z^ is the lidar altitude, K(A) is a calibration constant, 

2 

Q (,A,z,zl) is the two-way transmission between the lidar and height z, 
and D(z) can be any constant multiple of the atmospheric density 
profile. 
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In solving Eq. (2), the backscattered signal S(A,z) must first be 
obtained from the total detector output V(A,z) by subtracting the 
(internal plus external) background GCA), i.e., 

S(A,z) = V(A,z) - G(A ) . (3) 


Thereafter, the calibration constant K(A) is usually fixed by assuming 
that at the height z* where Eq, (2) attains its minimum value R^^CA ) , 
gas backscatterlng dominates particulate backscattering , implying that 
R m ^ n is close to unity (and hence can be estimated with relatively small 
error — see Section IV-D-4). 


This calibration procedure yields the expression from which R(A, z) 
is obtained in practice: 


R(A,z) = 


2 2 * * 
(z - z L ) S(A,z)Q (A,z L ,z )D(z ) 

* 2 * 2 
(z - z L ) S ( A. , z )Q (A,z l ,z)D(z) 


Rmin^ > 


(4) 


or 


where 


* * 
h(z - zl,z - zx J )s(A,z,z ) 

R(A,z) - — — ^min^^ > 

q(A,z,z )d(z ,z ) 


(5) 


* < z “ Z L> 

h(z - z L ,z - z L ) = — 

(z - z L ) 


2 * 


( 6 ) 


* S(A,z) 
s(A ,z,z ) = — 


* » 
S(A,z ) 


(7) 
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2 1 
A Q ( A , z l , z ) 


“ * “! 

Z 

r 

) 

q(A ,z,z ) = — — = exp ' 

±2 

I E(X,z')dz' 

| * 

Q » Z L» Z ) j 

i 

J 

z 

1 (8) 


* D(z) 

d(A )Z,z ) = — , 

D(z ) (9) 

and E(A jZ ) is the atmospheric extinction coefficient. In Eq. (8), the + 
(-) sign applies to an upward- (downward-) viewing lidar. 

Once R(A,z) has been obtained from Eq. (5), and Bg(X,z) has been 
obtained as described below, we can derive particulate backscattering by 
rearranging Eq. (1) to yield 

B p (A,z) = B g (.X,z)[R(X,z) - 1] . (10) 


In turn, the aerosol and cloud quantities required in climatic and 
environmental studies (e.g., particulate optical thickness, albedo, 
heating rate, mass) can be derived from B p (A,z) by using optical model 
conversion factors. As is well known, these conversion factors can vary 
widely when the model size, composition, and shape distributions are 
undefined; however, in some cases knowledge of B p (A,z) at two or more 
wavelengths might usefully constrain the possible range of optical 
models. 

In solving Eqs. (3) through (10), D(z) — and hence B ff (z)— is usually 

s 2 

obtained from nearby meteorological data; E(A,z) — hence Q (A,z,zl) — is 
usually provided by a model (often Improved by using lidar data); 

Rjiiin(^) is usually assumed to be unity or is estimated from previous 
aerosol data and other factors. Errors in these quantities and assump- 
tions, plus errors in the detector outputs V(X,z) and G(X) — hence in 
S(A,z) — produce corresponding errors in R(A,z), B p (A), and the other 
aerosol and cloud parameters derived from these results. Conventional 
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error-propagation analysis has shown that the resulting l<r relative 
uncertainties in R and Bp are given by 




where p and p are the absolute gas densities at heights z and z , 
respectively, 6x is the l<r uncertainty in each variable x. 


* * 
p = p(z ) 


(13) 


9 * 2 

C D d* is the covariance between D and D , and similarly for Cpp • 

[These covariances can usually be assumed to vanish except for |z - z* j 

<A Z » where A z is the vertical resolution of the input density profile. 

If z = z , then Cpj)* z = (<$D ) . However, if the density profile p(z) 

2 

contains a systematic bias error, then c pp* is nonzero at all heights. 
See Ref. 7 for details.] 

These results show that errors in determining the backscattered 
lidar signal [symbolized by 6s in Eqs. (11) and (12)] are but one of 
four sources of error in determining the retrieved cloud and aerosol 
quantities, R and B p . Errors in transmission (Sq), gas density (<$D), 
and lidar calibration (cSR^^) also occur and often dominate signal 
errors. Also, when B p /B g « R, uncertainties in retrieved B p , from any 
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of the four sources, are proportional to B g /B p . Thus, if other factors 
are equal, any increase of B p /B g (and hence of R) up to the point where 
B p /Bg » R tends to reduce errors in retrieved B p . 
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Ill SIMULATION PROCEDURE 


To test the uncertainty expressions of Eqs. (11) and (12), and also 
to develop insights and data analysis software, we developed a procedure 
(and computer code) to simulate the lidar measurement and retrieval pro- 
cess. The simulation procedure, described in detail elsewhere, ^ is 
illustrated in Figure 1. Given a set of inputs specifying lidar parame- 
ters, atmospheric backscattering and extinction profiles (gaseous and 
particulate), and background lighting, the code computes profiles of 
expected return signals and backgrounds. Then, random number generators 
inject random errors of appropriate size at each step of the measurement 
and retrieval sequence. The resulting profiles are compared to the 
input model profiles to see whether differences between the two are con- 
sistent with the algebraic uncertainty expressions, Eqs. (11) and (12). 
If desired, each simulation can be repeated an arbitrary number of 
times, using different random errors (drawn from the same distribu- 
tions); thereafter, root-mean-square (rms) differences between the set 
of retrieved profiles and the model profile are computed. The rms 
difference profile provides a more stable result for comparison with the 
algebraic lo- uncertainties, Eqs. (11) and (12). 

The simulation code also separates the algebraic uncertainties by 
source, as an aid in identifying the dominant sources for each mea- 
surement height and associated set of conditions. 
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FIGURE 1 SIMULATION PROCEDURE FOR EVALUATING LIDAR MEASUREMENT 
AND RETRIEVAL ERRORS 

Circles symbolize random number generators that inject simulated errors into derived 
quantities at appropriate steps. Numbers in parentheses indicate equations and tables 
in the text. 
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IV SIMULATION INPUTS 


A, Atmospheric Models 

The dependence of measurement and retrieval error on lidar 
wavelength, background lighting, density errors, and other factors is 
strongly determined by the particle and gas models used for input* For 
example, the wavelength dependence of particulate backscattering (along 
with other factors) determines whether it is advanageous to use alter- 
nate lidar wavelengths. Also, a strongly absorbing variable gas can 
produce significant transmission errors at its absorbing wavelengths. 
Moreover, model particle concentrations that give unrealistically large 
or small values for B p /Bg can yield unrealistically small or large rela- 
tive errors in retrieved B p [cf. Eq. (12)3. 

For these reasons we developed a set of input atmospheric models, 
each based on optical and physical measurements, as well as formation, 
transport, and removal processes that have been previously documented in 
the literature. In addition, we compared the resulting model back- 
scatter profiles at a wavelength of 0.694 pm with previous ruby lidar 
measurements, and ascertained that the model results were within the 
typical measurement range for the atmospheric regions and conditions 
they were intended to represent. For flexibility in simulating special 
situations, each model component can be multiplied by an arbitrary fac- 
tor over an arbitrary height range. 

The resulting atmospheric models have cloud and aerosol 
backscatter-to-extinction ratios, as well as wavelength dependences, 
that vary with height, latitude, and particle source (e.g., volcanic, 
Saharan, marine) in accordance with the assumed particle size distribu- 
tion, composition, and shape. The models, which include molecular 
density, ozone, cirrus and noctilucent clouds, and mesospheric, stratos- 
pheric, upper tropospheric, Saharan, and marine aerosols, are described 
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in Appendix A, Figures 2(a) through 4(a) show the composite model pro- 
files of backscatter mixing ratio (Bp/Bg) that we will use for the simu- 
lations in this paper. The associated signal profiles are discussed in 
Section V-A. 

As explained more fully in Appendix A, the inclusion of mesospheric 
aerosol layers in our models is not intended to suggest that such layers 
are always, or even typically, present. Rather, the intent is to con- 
struct model layers that are consistent with several previous mea- 
surements (and transport processes), and then to see how well the orbit- 
ing lidar could detect such layers if indeed they were present. Also, 
we have emphasized marine aerosol models over continental ones because a 
major strength of any orbiting lidar will be its ability to fill data 
gaps over the world's vast ocean areas. Note that the marine aerosol 
models do include particles of continental origin (see Appendix A). 


B. Lidar Parameters 

The intent of this study is to simulate measurements made by a 
first-generation spacebome lidar, flying on the space shuttle and using 
straightforward measurement techniques with a history of analogous mea- 
surements on earth. Therefore, we adopted the parameters of the base- 

4 

line shuttle lidar determined in a recent facility definition study by 
General Electric. These parameters, listed in Table 1, satisfy the fol- 
lowing criteria: 

• They are compatible with the weight, power, safety, and other 
restrictions for shuttle science payloads. 

• They satisfy a -set of preliminary eye safety guidelines 
developed by General Electric.* 


These guidelines assume the ANSI criteria for maximum permissible 
occular exposure, viewing with the naked eye by day and with a 10-inch 
telescope by night, a safety factor of 10 for laser beam scintilation 
effects, and an additional factor of 2 or 3 to account for the differ- 
ence between average and peak laser energies in Gaussian or multimode 
laser beams, respectively. 
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Table 1 


ASSUMED LIDAR PARAMETERS 



Wavelength, X (pm) 

Parameter 

1.064 

0.532 

0.355 

Transmitter 




Pulse energy (J, simultaneous) 

1.1 

0.55 

0.20 

Repetition rate (s *) 

10 

10 

10 

Full beamwidth (mr) 

5*. 0.V 

5*, 0.4 1 * 

•k ■+■ 

5 , 0.4 T 

Receiver 




Diameter (m) 

1.25 

1.25 

1.25 

Full beamwidth (mr) 

5.5*, 0.5 1 * 

5.5*, 0.5 1 * 

5.5*, 0.5 f 

Filter bandwidth (nm) 

4.0*, 0.4 1- 

0.2*, 0.02 t 

4.0*, 0 .4 1 * 

Filter transmission 

0.50*, 0.25 + 

0.50*, 0.25 + 

0.20*, 0.10 + 

Overall transmission 

0.34*, 0.17 1, 

0.25*, 0.12 1 * 

o.io*, 0 . 05 * 

Quantum efficiency 

0.04 

0.20 

0.30 

Detector dark counts (ys"^) 

0.03 

0.03 

0.03 

Detector gain 

10 6 

10 6 

10 6 

Altitude, 300 km 





Nighttime configuration. 
^ Daytime configuration. 
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• All three wavelengths, though transmitted simultaneously, are 
obtained from a single Nd:YAG laser i This type of laser has 
been highly engineered recently for flight applications, rugged- 
ness, and efficiency, 

• No special tuning, line narrowing, or wavelength stability is 
required for either the transmitter or the receiver. 

In stam, the lidar parameters define a baseline spaceflight system. 
Nevertheless, certain engineering issues remain to be resolved regarding 
frequency tripling and automatic changing of receiver field of view and 
filter bandwidth when crossing from dark to daylit areas and back. 


C. Background Lighting 


To give an idea of the wide range of possible background lighting 
conditions, we have chosen three sets of scenarios: sunlit cloud, sun- 

lit ocean, and full moonlit cloud, each with the illuminating source 
(sun or moon) at zenith angle 0 o . For our purposes both clouds and the 
ocean can be modeled as Lambertian reflectors — i.e., reflected radiance 
I B (X) can be modeled as independent of the angle of emergence. This 
yields, for the spectral radiance I B (X) emerging toward the zenith, 


2tt 


it/2 


I B (*) = A 


F(A )cos0 o // a f 


cos© sin0 


= A F( A ) cos0 o /(n sr) 


(14) 


where A is the cloud or ocean albedo and F(A ) is the solar (or lunar) 
spectral irradiance Note that we do not decrease I B (X) by a clear- 
atmosphere, two-way transmission or increase Ig(A) by the clear- 
atmosphere reflectivity. These two factors tend to cancel; moreover, 
their inclusion is not warranted, given the approximations already made 
in the value and spectral and angular dependence of cloud and ocean 
reflectivity. (Also, our simulation results below 20 km are not 
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sensitive to changes in background by factors of 3 or less.) Table 2 
shows results of Eq. (14) for the three chosen scenarios. 

Table 2 

BACKGROUND LIGHTING (ZENITH UPWARD SPECTRAL RADIANCE) 

FOR A DOWNWARD-VIEWING LIDAR AND SEVERAL SCENARIOS* 



I B 

(mW m ^ sr” ^ nm'^) 

X 

Sunlit Cloud 

Sunlit Ocean 

Moonlit Cloud 

(pm) 

(A = 0.8) 

(A = 0.1) 

(A =» 0.8) 

1.064 

168 

22 

168 x 10" 6 

0.532 

460 

58 

460 x 10~ 6 

0.355 

280 

34 

280 x 10“ 6 


Calculated from Eq. (14) assuming the solar spectral irradiances 
F(x) of Thekakera^ and cos 0=1. 


D. Error Sources 
1 . Signal 

In most cases the wide dynamic range of backscattered signal pro- 
files received by a spaceborne lidar will produce a detector output that 
must be measured partly in the pulse-counting mode and partly in the 
current mode; the resulting two profile parts must be accurately spliced 
together (see Section V-A-l , V-B-l ) . In the pulse- counting mode the lo- 
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uncertainty in the number of signal counts, N g (per range element), is 
given by 



«N S V N S + N B + N D 

N s N S 


(15) 


where Ng and Ng are the number of output counts arising from background 
lighting and internal detector noise, respectively. [This assumes that 
N B + Ng is determined accurately enough that its uncertainty is negligi- 
ble compared to Eq. (15). See Ref. 7 for details.] In the current mode, 
errors are caused by factors such as detector shot noise and saturation, 
amplifier noise and nonlinearity, and digitizer truncation and satura- 
tion. If gain switching or logarithmic amplifiers are used to compress 
the dynamic range of detector output, additional error sources are 
Introduced. Finally, matching the counting and current parts of the 
output profile becomes another error source. These errors must be 
reduced to acceptable levels if detector outputs from the spaceborne 
lidar are to be useful over all height regions of interest. This prob- 
lem has been solved for certain height ranges by using a variety of 
techniques in many ground-based and airborne lidar systems. However, 
each technique has its own error characteristics, and it is not clear 
which technique, or possibly a new technique, ^ would be used for a 
shuttle lidar system. 

In our simulation program we have confronted this simulation prob- 
lem in the following manner. At each measurement height, N s is calcu- 
lated from the lidar equation as shown in Ref. 7, and Ng and Ng are cal- 
culated from the background-lighting, receiver, and detector charac- 
teristics. The pulse-counting value for 6 S/S, Eq. (15), is then calcu- 
lated. As Ng increases, this result decreases toward zero. However, it 
is clear that increasing Ng also leads to pulse overlap, making pulse- 
counting impossible, and forcing signal measurement into the current 
mode. At that point, the current-measurement and transition errors men- 
tioned above take over, and Eq. (15) is no longer a valid expression of 
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total signal uncertainty. Since no universally valid expression 
describes the current-mode errors mentioned above, our program uses 


<$S 

= MAX 

S 




( 16 ) 


That is, whenever the photon-counting expression becomes less than a 
prespecified minimum, (^S/S) m ^ n , that minimum is used for (6S/S). This 
assumes that the shuttle lidar hardware can measure the current-mode 
part of the output profile and splice it to the pulse-counting part with 
the accuracy (<$S/S) m £ n . To highlight the importance of this assumption, 
the program plots profiles of expected detector output (backscattered 
signal and background) as a function of height, in terms of both pulse 
rate and annode current, relative to certain detector benchmarks, such 
as the point where pulse overlap becomes highly likely and the point 
where current measurements become sufficiently accurate. (The bench- 
marks are the 1- and 100-count ys ^ lines in Figures 2 through 4.) 

These plots can then be used in designing signal-processing and record- 
ing hardware that satisfy the (cSS/S)^.^ requirement. 

For the simulations reported here, we used 



Although this requirement is stringent, we feel it is attainable using 

9 9 

the technique described by Evans. Relaxing this requirement to 0.01 
would have little effect on the aerosol and cloud simulations shown in 
this part; however, it would significantly increase the minimum density 
and temperature errors shown in Part 3. (See also Section V-A of this 
part.) 
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Finally, we require an expression for <5s/s to use in Eqs. (11) and 
(12). From Eq. (7) it follows that 


where 

£ it 

S = S(A, Z ) . (19) 

2. Molecular Density 

The problem of obtaining molecular density information from 
conventional meteorological sources (radiosondes, satellites, and rock- 
etsondes) is discussed in Ref. 7. From a radiosonde ascent, density 
profiles can be determined with an accuracy of about 1 percent. How- 
ever, since most shuttle lidar data will come from remote regions of the 
globe (e.g., the large ocean areas not covered by radiosondes), one must 
assume that current radiosonde density data will not be available for 
most shuttle lidar data. The alternative sources of conventional infor- 
mation are: 

• Carefully selected seasonal and regional model density profiles. 

• Density profiles obtained from current satellite and radiosonde 
network analyses and historical rocketsonde data using hydro- 
dynamically consistent interpolation methods. 

The quality of such density profiles has recently been evaluated by 
National Weather Service personnel in connection with the SAM II and 
SAGE satellite programs (which require density profiles in remote areas 
for input to routine data analysis). These evaluations (Ref. 23, Appen- 
dix D, by J. Laver) show that expected density errors depend in a fairly 
complicated way on height, the time difference between the lidar mea- 
surement and the conventional analyses, the hemisphere of the lidar mea- 
surement, and even the latitude within the southern hemisphere. How- 
ever, for purposes of this study, representative values can be 
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summarized as in Table 3* The indicated 5D/D values were used for the 
simulations reported here. The expression given by Ref. 7 was used for 

C DD*« 

Table 3 

REPRESENTATIVE lcr UNCERTAINTIES IN DENSITY PROFILES 
DERIVED FROM CONVENTIONAL SOURCES 
Values apply to densities defined at 
geometric (not pressure or geopotential) height 


Data Source/z 

0-30 km 
(%) 

30-50 km 
(%) 

50-90 km 
(%) 

Radiosonde measurement 
within 100 km and 6 h 
(no intervening frontal 
activity) 

1 



Hydrodynamically consistent 
interpolation of current 
radiosonde and satellite 
data, plus historical 
rocketsonde data, in data- 
sparse area 

2 

3 to 5 

10 

Regional and seasonal 
models 

3 

5 

10 


The uncertainties used for the simulations in this paper are: 

Low-latitude : 2% below 30 km, 3% above. 

Mid-latitude : 3% below 40 km, 10% above. 

High-latitude : 3% below 40 km, 10% above. 

3. Transmission 

As shown in Ref. 7, transmission uncertainties for the wavelengths 
of interest are given by 
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* 

where t x (A,z,z ) is the optical thickness of constituent /mechanism x 
between heights z and z*; and p, sg, and 3 denote (cloud and aerosol) 
particles, total gaseous scattering, and ozone absorption, respectively. 
(As recently pointed out by Young, x__. includes Raman and Brillouin 
scattering as well as Rayleigh scattering.) For the reasons cited in 
Ref. 7, we use 

6T3(A,z,z ) = 0.2t3(A,z,z ) • (21) 


Ref. 7 uses the conservative overestimate 

* * 

5 Xg g(AjZ,z ) 0. lx gg(A , z ,z ) 


( 22 ) 


for gaseous scattering optical thickness uncertainties. However, this 
overestimate becomes important at the short wavelength, 0.355 ym. 
Therefore, we now use the more exact expression 


^ r ^ ^ *] ^ 

'5 1 't sg (*,z,z ) = [SAPCz.z )/AP(z.z )J T sg( x > z > z ) 


(23) 


where ^P(z,z ) is the pressure difference between z and z . In lidar 

data analysis the pressure difference is obtained from the molecular 

density profile, and thus the pressure difference uncertainty, 

* 

^Ap(z» 2 ) , is a function of the molecular density uncertainty profile, 

it 

6p(z). Our simulation program obtains SAP(z>z ) from the p(z) and <5p(z) 
profiles using error-propagation techniques described in the companion 
paper 


For the reasons described in Appendix B, we retain the expression 

1/2 


8x 0 (A,z,z ) = 0.5 


2 2 
£ T ap( ^» Z » Z )J ”*■ ^ cp ( ^ y z I z )j 


( 24 ) 


used in Ref. 7. 
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Note that Eq. (20) includes no term to describe NO 2 absorption 
effects. We intend to add such a term in the future and to include typ- 
ical NO 2 profiles in our atmospheric models as better data become avail- 
able. However, for the present we have calculated the error incurred by 

neglect of NO 2 in Eq. (20). As inputs, we used the NO 2 profiles summar- 

25 26 

ized by Noxon^ and the absorption coefficients reviewed by Hudson. 

Taking a worst case of nighttime midlatitude conditions, with an uncer- 
tainty in NO 2 concentrations of ±100 percent between 5 and 20 km and ±50 
percent between 20 and 35 km, led to the following results for (<$q/q)^Q 2 
at wavelength 0.355 pm: 


z 

(km) 

* 

z 

(km) 

(Sq/q) 

no 2 

5 

8 

0.0002 

20 

8 

0.0008 

28 

8 

0.0016 

35 

8 

0.0026 


These values are small compared to the other terms in Eqs. (20), (11), 
and (12). Moreover, for daytime, lower latitudes, and the longer 
wavelengths used in this study, (5q/q) N0 2 values are even less.^ 


^min 

Reference 7 discusses the problem of determining expectation values 
<R m i n > and uncertainties as a function of latitude, lidar analysis 

height range, volcanic effects, and other factors. That study shows 
that <R tnin > and < ^R tn ^ ri can be estimated from the large set of particle 
number profiles measured by balloonborne particle counters at sites dis- 
tributed from 90°S to 84°N in the past 12 years, provided the particle 
number data are combined with realistic models of relative size distri- 
bution and composition. As an example, Ref. 7 presents values of 
and 5R using model size distributions and compositions that are 
appropriate for particles in and near the peak of the stratospheric 
aerosol layer. 
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Since that work was done, it has become clear that significantly 
different types of size distribution, and different compositions, should 
be used to describe aerosols in the upper troposphere, where R m ^ n is 
frequently found* Using such a height-dependent optical model with the 
particle number profiles tends to increase the results for < R m -£ n > and 
6R to i n found in Ref. 7. While the scope of this study does not permit us 
to derive a complete new set of values to replace those in 

Ref. 7, we have nonetheless performed sample calculations suggesting 
that the revised optical model increases the < R m £ n “ 1> and ^'r^ values 
of Ref. 7 by about a factor of 2. (However, in northern high latitudes, 
the values in Ref. 7 practically unaffected, because R m ^ n rarely occurs 
in the troposphere for such latitudes.) For the current simulation 
study, we feel that it is sufficiently accurate to use - 1> and 

^Rjnin values that are double those of Ref. 7, with the high-latitude 
exception just noted. However, because of the importance of R min to all 
lidar results derived by the scattering ratio/normalization method [cf. 
Eq. (5)] , we emphasize that analyses of actual lidar data should take 
care to estimate the best possible values for and for the 

conditions under which the data were taken. 

The results in this section reinforce a conclusion of Ref. 7, 
namely that a carefully estimated value for < R m -^ n > is in many cases sig- 
nificantly more accurate than the common assumption of R m ^ n = 1, and 
can significantly improve the accuracy of retrieved lidar data products 
over large height ranges. In fact the results of this section suggest 
that the error in using R m ^ n = 1 is about double that previously 
thought. 
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V SIMULATION RESULTS 


A. Signal and Background Profiles 

Figures 2 through 4 show single-shot profiles of signal and signal 
plus background for the lidar parameters of Table 1 and the three asso- 
ciated model atmospheres each for a nighttime and a daytime condition. 
The detector outputs have been expressed both as photoelectron counts 
per microsecond and as microamperes, assuming a detector internal gain 
of 10^, which is typical for a commercially available photomultiplier 
recently tested for lidar applications by Evans. The solid vertical 
line in each plot indicates the counting rate at which output pulses 
(each ~20 ns long ) become likely to overlap, thus causing pulse- 
counting errors. The dashed vertical line indicates the minimum output 
for which conventional current measurements can be made with the accu- 
racy specified by Eq. (17). 

In all cases the dynamic range of the total detector output (signal 
plus background) in the current-measurement range is less than three 
orders of magnitude (although the return from the surface or opaque 
clouds would be considerably larger than the atmospheric returns shown). 
This behavior contrasts with many ground-based or airborne lidars, which 
can encounter many decades of signal decrease from the range-squared 
effect alone. In contrast, for a shuttle lidar altitude of 300 km, 
range-squared effects vary by less than a factor of 3 between 100 km and 
the surface, and actually counteract the increase of atmospheric 
density. 

The major problem in accurate measurement of the output profiles is 
that, at night, all profiles span the transition region between pulse- 
counting and current measurement. Because of the normalization method 
assumed here to calibrate received signals In terms of backscattering 
coefficient, the data acquisition method must accurately splice together 
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the pulse and current portions of the profile in this transition region. 
A method for doing this was recently described by Evans, ^ b ut to our 
knowledge it has not been implemented. As noted in Section IV-D-1 , we 
have assumed that current measurements and a pulse-to-current splice can 
be made with an accuracy of 0.5 percent. Although we feel that this 
goal can be met, it is worth emphasizing that this assumption implies a 
certain amount of hardware development. 


B. Retrievals 

1 • Tropical , Nonvolcanic , Cloud-Free , Saharan 

a. Nighttime ( Zenith Moonlit Cloud ) , Ax » 200 km 

The model atmosphere used in this scenario [Figure 2(a)] includes 
no clouds; however, we have used a background radiance appropriate for 
clouds illuminated by a full, zenith moon (Table 2) to make our results 
conservative. As can be seen from Figure 3(a), 1.06-pm signal-to-noise 
ratio [Eq. (15)] above 20 km would improve by assuming less background 
(e.g., a moonlit ocean); however, the 0.53-jjm signal-to-noise ratio is 
signal-dominated for all heights below 35 km, and hence is insensitive 
to this assumption. 

Figure 5(a,b) shows simulated retrievals for backscatter mixing 
ratio, assuming the error sources given by Eqs. (15)-(18), (20)-(24), 
Table 3, and Section IV-D-4, with shots integrated for 28 s (total 280 
shots), yielding a horizontal resolution of -200 km for a spacecraft 
velocity of -7 km s * . Vertical resolution, shown by the spacing of the 
dots, is 0.5 km below 20 km and 1 km above, except above -29 km where 
range bins have been combined to improve the signal-to-noise ratio. The 
solid curves give the model (i.e., "exact") profile; the dots give the 
lidar retrieval; and the error bars give the expected Icr uncertainty 
calculated from Eqs. (11) or (12). Although for this case the error 
bars somewhat overestimate the rms difference between the retrieval and 
the model, close inspection has shown that this was caused by fortui- 
tously small random errors drawn for R min and transmission in this case. 
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FIGURE 5 LOW-LATITUDE NIGHTTIME SIMULATION RESULTS, 1.064 AND 0.532 jum, 
USING CONVENTIONAL DENSITY DATA 

(a,b) Backscatter mixing ratio profiles; (c,d) Relative uncertainty in particulate 
backscattering broken down by source. 
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Many simulations run in this study and previously^ have substantiated 
the validity of the algebraic expressions Eqs. (11) and (12), 

For this scenario both wavelengths give quite good retrievals below 
7 km and between 22 and 30 km. Between 8 and 20 km, where backscatter 
mixing ratios are near their minimum, relative uncertainties increase, 
and the quality of the retrievals is poorer, especially for A = 0.53 ym. 
The reason for this is shown in Figure 5(c,d), where the relative uncer- 
tainty in particulate backscattering is broken down by source. For 0.53 
ym, at all heights above 6 km, uncertainties in molecular density are 
the dominant source of error. This simulation assumes a relative den- 
sity uncertainty of 0.02 below 30 km, which is appropriate for density 
profiles in remote areas derived by interpolation from current conven- 
tional data. Because the backscatter mixing ratio B^/B^ at 0.53 ym is 
only about 0.01 to 0.02 between 8 and 20 km, the density uncertainties 
alone lead to B^ uncertainties of 100 to 200 percent in this region. At 
1.06 ym, backscatter mixing ratios are considerably larger (~0.08 to 
0.20), and the density-induced relative uncertainty in B^ is propor- 
tionately reduced. 

Note, however, the greatly increased calibration uncertainty at the 
longer wavelength, caused by uncertainties in W Uncertainties in 
Rni n (1.06 ym) considerably exceed those in ym) , because even 

very tenuous aerosol layers can perturb ym) by several per- 

cent, whereas the same aerosol layers perturb ^^(0-53 ym) much less. 
(Compare the backscatter mixing ratios in Figure 2.) Thus, at 0.53 ym, 
one is more likely to find a "calibration layer" where gas backscatter- 
ing strongly dominates particle backscattering, thus permitting one to 
calibrate the entire received signal profile in terms of backscattering 
coefficient. 

Between 7 and 5 km, in the top of the strong Saharan dust layer, 
relative uncertainties are greatly reduced at both wavelengths because 
of the sharp increase in B^. However, with increasing penetration into 
the Sahara layers, transmission uncertainties begin to degrade the 
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retrieval, and the overall uncertainty again increases. These transmis- 
sion errors are discussed in more detail in connection with the next 
scenario • 

Note that signal errors do not make a significant contribution in 
this scenario, except above -30 km, especially at 1,06 ym. However, the 
weaker signals at these heights have had the indirect effect of forcing 
a coarser vertical resolution there (through rebinning to reduce signal 
errors). This effect would have become more serious, and penetrated 
lower, if the integration time (horizontal resolution) had been 
decreased below 28 s (-200 km). 


b • Daytime ( Zenith Sunlit Ocean ) , A x » 100 km 

Figure 6 shows simulated retrievals and uncertainty breakdowns for 
the same model atmosphere as Figures 5 and 2(a), but with the lidar in 
its daytime configuration (Table 1) and a background radiance that 
approximates that of a sunlit ocean (Table 2). (The sun is assumed to 
be near the zenith, but not exactly at the zenith, as this would cause a 
very large specular reflection into the lidar telescope.) As can be 
seen from Figure 3(b), signals above 20 km are background-dominated, and 
no useful retrievals are possible; hence, results are shown only below 
20 km on an expanded scale. Note, however, that within the Sahara layer 
and the underlying marine layer, even the daytime background is small 
compared to the backseat tered signal at both wavelengths. Hence, both 
vertical and horizontal resolution can be fine. This simulation assumes 
Ax « 100 km (140 shots, 14 s), with Az = 0.25 km below 7 km, and Az = 1 
km above. (Note, however, the coarser vertical resolution above 10 km 
for 1.06 ym, caused by rebinning.) The only reason for making A?^ as 
large as 100 km is to provide a relatively accurate signal profile in 
the clean upper troposphere, for normalization (calibration) of the 
entire profile. Measurement of the vertical structure of the Saharan 
and marine aerosol layers, and of their relative backscattering 
strengths (to within the transmission errors noted below) is possible on 
a single-shot basis. This also applies at night. 
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6Bp/Bp — percent 

FIGURE 6 LOW-LATITUDE DAYTIME SIMULATION RESULTS, 1.064 AND 0.532 pm, 
USING CONVENTIONAL DENSITY DATA 

(a,b) Backscatter mixing ratio profiles; (c,d) Relative uncertainty in particulate 
backscattering broken down by source. 


54 





Figure 6 shows that retrievals at both wavelengths are very good in 
the upper part of the Saharan layer, but that retrievals are degraded by 
transmission errors as penetration increases. However, because 
transmission errors change the shape of retrieved profiles in a gradual 
way, the presence and height of the lower Sahara layer at 2 km and the 
marine layer near the surface are well detected, even though backscatter 
strengths are not retrieved correctly. Thus, for example, an inference 
of marine aerosol layer depth should be possible in spite of transmis- 
sion errors. This inference holds as long as extinction is not large 
enough to degrade signal measurement accuracy significantly. 

A final point to notice in Figure 6 is that the error bars in Fig- 
ure 6(a), in contrast to those in Figure 5(a), appear to underestimate 
systematically the lo- difference between the retrieval (dots) and the 
model (solid line). Also, the retrieved profile is systematically 
offset to the right (i.e., too large), even above the height where 
transmission errors systematically distort the profile. The reason for 
the offset is a relatively large (_20 percent) negative error in the 
signal profile at 9.25 km, caused by drawing a random number that dif- 
fered from the mean by about twice the distribution standard deviation, 
SS = 0.10S . The normalization procedure [Eq. (4)] then led to the 
systematic offset of the entire profile by an amount that exceeded the 
1 er uncertainty 6S(z ). This type of error will also occur occasionally 
in analyzing actual lidar data, and will then produce retrieved profiles 
that differ systeraa tically from the true profile by more than the 
derived lo- error bars. Note that the lcr uncertainties given by Eqs. 

(11) and (12) describe the expected rms deviation of retrievals from 
actuality in an ensemble containing many profile measurements, rather 
than within one retrieved profile. 
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2 


Midlatitude, Volcanic , Cirrds , Marine 


a . Nighttime ( Zenith Moonlit Cloud ) , /\x « 200 km 

Figure 7 shows simulated retrievals for a midlatitude model atmo- 
sphere that includes a mesospheric aerosol layer at 46 to 54 km (see 
Appendix A). This simulation assumes density uncertainties of 3 percent 
below 40 km and 10 percent above. Even the relatively large density 
uncertainty of 10 percent is small compared to the mesospheric aerosol 
backscatter mixing ratio peak [~ 1.0 at 1.06 ym; ~ 0.6 at 0.53 ym; see 
also Figure 2(b)]. Thus, density uncertainties do not prevent a mea- 
surement of the mesospheric aerosol at either wavelength. However, sig- 
nal errors do prevent a useful measurement above 45 km at 1.06 ym; 
hence, the mesospheric layer is not detectable at this long wavelength. 
At 0.53 ym, signal errors are considerably less, because of both 
increased backscatter coefficient and increased quantum efficiency. As 
can be seen in Figure 7(b,d), the peak of the mesospheric aerosol layer 
is unambiguously detected at this wavelength, with a measurement uncer- 
tainty of about 25 to 50 percent for the aerosol model used here (see 
Appendix A. 3). 

In the stratospheric aerosol layer (16 to 25 km) the moderate vol- 
canic enhancement modeled here yields backscatter mixing ratios that are 

a factor of 3 to 6 larger than those of the background midlatitude stra- 

16 18 

tospheric aerosol. » As a result, the relative retrieval error caused 
by density errors is proportionately reduced. Both wavelengths give 
good retrievals in the stratospheric aerosol peak; however, the region 
of good retrieval extends farther from the peak at 1.06 ym because of 
the reduced effect of density errors at that wavelength. 

Both wavelengths give good quantitative retrievals of the subvisi- 
ble and visible cloud layers, and both detect the presence of the 
elevated and surface marine aerosol layers below the clouds. However, 
transmission errors degrade the quantitative retrievals below the lowest 
cloud . 
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FIGURE 7 MID-LATITUDE NIGHTTIME SIMULATION RESULTS, 1.064 AND 0.532 
USING CONVENTIONAL DENSITY DATA 

(a,b) Backscatter mixing ratio profiles; (c,d) Relative uncertainty in particulate 
backscattering broken down by source. 
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b 


Daytime ( Zenith Sunlit Cloud ) , « 100 km 


Figure 8 shows results for the same model atmosphere as Figures 7 
and 2(b), but now for the daytime lidar configuration (Table 1) and a 
sunlit cloud background radiance. At both wavelengths, signal errors 
prevent useful measurements above 35 km, so the mesospheric aerosol 
layer is not detectable during the daytime. However, at 0.53 ym, signal 
errors are small enough that the moderate volcanic stratospheric aerosol 
layer is barely detectable; notice, though, that use of noisy signals 
above 25 km for normalization leads to a normalization error that 
affects the uppermost cloud retrieval significantly. At 1.06 ym, even 
the moderate-volcanic stratospheric aerosol layer is not detectable by 
day, and normalization is achieved below its height, above and between 
the clouds. 

As with the nighttime case, both wavelengths give good quantitative 
retrievals for the backscattering coefficient at the top of each sub- 
visible or visible cloud layer, with the exception noted above caused by 
poor normalization. However, transmission errors in the lowest cloud 
significantly degrade quantitative retrievals below that cloud's top. 

As noted previously, these errors do not affect the heights inferred for 
the lowest aerosol layers, but they do prevent quantitative inferences 
of the amount of aerosol present. Notice that heights are retrieved 
with a vertical resolution of 0.25 km in this simulation; in fact 
signal-to-noise ratios permit even finer vertical resolution — see 
Section V-B-3-b. 

Notice that the transmission errors in Figure 8(a,b) are opposite 
in sign to those in Figure 7(a,b). This simulates the use by the 
analyst of two different erroneous transmission profiles in the two dif- 
ferent scenarios: overestimating the transmission in Figure 7(a,b), and 

underestimating it in Figure 8(a,b). 
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FIGURE 8 MID-LATITUDE DAYTIME SIMULATION RESULTS, 1.064 AND 0.532 jum, 
USING CONVENTIONAL DENSITY DATA 

(a,b) Backscatter mixing ratio profiles; {c,d) Relative uncertainty in particulate 
backscattering broken down by source. 
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3, High-Latitude , Nonvolcanic , Noctllucent Cloud , Marine 


a . Nighttime ( Moonlit Cloud , cos 0 Q = 10 . 5^) , *• 2000 km 

Figure 9 shows results for a high-latitude model atmosphere [Figure 
2(c)] that includes both a noctilucent cloud (at 80 km) and an enhanced 
mesospheric aerosol layer (63 to 73 km) of the type that accompanied 
noctilucent clouds (NLCs) in the measurements of Fiocco and Grams. 

This simulation assumes density uncertainties of 3 percent below 40 km 
and 10 percent above. As in the previous case, density errors are not 
critical for measurements of the mesospheric layers, but signal errors 
are, especially at 1.06 ym. To see whether increased measurement time 
might significantly aid the 1.06-ym measurement, this scenario assumes 
integration for 280 s (2800 shots, ^x w 2000 km). At the high- latitude 
extreme of an orbit with inclinations 70 to 80°, such a distance on the 
suborbital path would cover only _5° of latitude and 40° of longitude; 
thus, the mesospheric aerosol layer (though probably not the NLC) may be 
uniform enough in this integration distance to provide a meaningful 
average • 

Figure 9(a) shows that even with such a long integration period 
1.06-ym signal errors are too large to permit a useful retrieval of the 
mesospheric aerosol layer. However, because the 1.06-ym signal from the 
noctilucent cloud is comparable to that from the stratospheric aerosol 
[see Figure 3(c) 1, the NLC is readily detected. Figure 9(b) shows a 
simulated retrieval at 0.53 ym for the same integration period. The 
retrieval of both the mesospheric aerosol and the NLC is quite good, and 
Figure 9(d) shows that, for this wavelength and integration time, signal 
errors are a minor source of uncertainty. In fact, a similar retrieval 
could be obtained with an integration time of 70 s (700 shots, ^z » 500 
km). At the 0.53 ym wavelength, the NLC layer can be detected with a' 
single shot. [See Figure 3(c)]. Hence, the layer's height could be 
monitored on a very fine horizontal scale, for example to trace out wave 
motions. (At 10 Hz and 7 km s , the suborbital shot-to-shot distance 
is 700 m.) 
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FIGURE 9 HIGH-LATITUDE NIGHTTIME SIMULATION RESULTS, 1.064 AND 0.532 /um, 
USING CONVENTIONAL DENSITY DATA 

(a,b) Backscatter mixing ratio profiles; (c,d) Relative uncertainty in particulate 
backscattering broken down by source. 
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In the stratosphere, the high-latitude nonvolcanic aerosol layer is 
lower (~12 to 18 km) than at other latitudes (cf. Figure 2), and the 
increased gas backseat tering B at the lower heights yields smaller 

O 

backscatter mixing ratios B /B for a given particulate backscattering 

t o 

Bp. At 0.53 urn, the model stratospheric aerosol backscatter mixing 
ratios (0.03 to 0.05) are comparable to the assumed lo* uncertainty in 
molecular density (0.03). Hence, detection of the layer is barely pos- 
sible, and measurement errors make it impossible to retrieve the layer's 
vertical structure. At 1.06 ym, the considerably larger backscatter 
mixing ratios (0.2 to 0.3) greatly reduce the effect of density errors, 
permitting a good retrieval of layer vertical structure. Uncertainty in 
R min i ncreases the uncertainty of absolute backscatter retrievals to 30 
to 60 percent [Figure 9(c)], but this uncertainty does not affect 
retrieved vertical structure. 

Both wavelengths give good quantitative retrievals of the elevated 
marine aerosol layers (1 to 3 km). In fact, as the next scenario shows, 
considerably finer vertical resolution is possible, permitting a 
retrieval of the surface-based marine aerosol layer's height. 

b. Daytime ( Sunlit Cloud or Ice , cos 0 Q = 0,5 ) 

Ax ~ 100 km 

Figure 10 shows results for the same model atmosphere as Figures 9 
and 2(c), but with the daytime lidar configuration (Table 1) and a back- 
ground approximating a sunlit cloud (Table 2). At both wavelengths the 
nonvolcanic stratospheric aerosol is retrieved very poorly (if at all); 
however, the marine aerosol layers are retrieved very well with a verti- 
cal resolution of 0.1 km. 
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FIGURE 10 HIGH-LATITUDE DAYTIME SIMULATION RESULTS, 1.064 AND 0.532 jum, 
USING CONVENTIONAL DENSITY DATA 

(a,b) Backscatter mixing ratio profiles; (c,d) Relative uncertainty in particulate 
backscattering broken down by source. 


63 






VI SUMMARY AND CONCLUSIONS 


We have evaluated the expected accuracy of retrieved particle back- 
scatter information for a proposed shuttlebome lidar operating over 
three model atmospheres, each with nighttime and daytime background 
lighting. The model atmospheres, representing low, middle, and high 
latitudes, include profiles of total gas density, ozone, tenuous clouds 
(subvisible and visible), and mesospheric, stratospheric, and tropos- 
pheric (Saharan and marine) aerosols. Each model is based on previous 
optical and physical measurements, so as to represent conditions likely 
to occur in practice. The background lighting conditions simulate radi- 
ances from sunlit clouds, a sunlit ocean, and moonlit clouds. 

The data profile from each lidar wavelength was assumed to be 
analyzed independently of data from other wavelengths, using the 
scattering ratio/normalization method of analysis. A further assumption 
was that only conventional total gas density profiles were available, 
and that these had rms errors ranging from 2 to 3 percent below 30 km to 
5 to 10 percent above 40 km in data-sparse regions of the globe, where 
most shuttle lidar measurements will be made. For each scenario, 
expected errors in retrieved particle backscatter were calculated from 
algebraic expressions and also checked by numerical simulations using 
random number generators. Both the algebraic expressions and the simu- 
lations include all four sources of error that we have encountered in 
actual lidar measurements and analyses: namely, errors in (1) measuring 

the backscattered signal, (2) accounting for transmission losses, 

(3) determining the total gas density profile at the lidar location, and 

(4) calibrating signals in terms of backscatter using an atmospheric 
layer where gas backscatter dominates particle backscatter. 
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The results show that useful retrievals can be made for the follow- 
ing constituents and conditions: 

• By day : Vertical structure of tenuous clouds (subvisible and 

visible), Saharan aerosols, and boundary-layer aerosols (at both 
0.53- and 1.06-ym wavelengths), and strong volcanic stratos- 
pheric aerosols (at 0.53 ym). Quantitative backscatter can be 
determined for all of the preceding, provided particulate opti- 
cal depth does not exceed ~0.3. 

• night : Vertical structure and quantitative backscatter for 

all of the above, plus upper tropospheric and stratospheric 
aerosols (at 1.06 ym) and mesospheric aerosols and noctilucent 
clouds (at 0.53 ym). (Again, quantitative backscatter 
retrievals suffer when particulate optical depth exceeds -0.3.) 


These results were obtained with vertical resolutions of 0.1 to 0.5 
km in the troposphere and 0.5 to 2.0 km in the stratosphere and above, 
except within the mesospheric aerosol and noctilucent cloud layers, 
where resolutions of 1.0 and 0.25 km, respectively, were used. Horizon- 
tal resolution was 100 km in the troposphere, 200 km in the stratosphere 
and lower mesosphere (below -60 km), and 2000 km in the upper mesosphere 
(to -85 km). Vertical structure of the stronger layers (noctilucent and 
cirrus clouds, lower tropospheric aerosols) could be observed on much 
finer horizontal scales (sometimes on single shots — about 700 m apart); 
moreover, these fine-scale results could be converted to quantitative 
backscatter if a transmit energy monitor or stable transmit energy were 
provided . 


In general, the principal virtue of the 1.06-ym wavelength is its 
reduced sensitivity to gas density errors, which are the dominant error 
source for upper tropospheric and nonvolcanic stratospheric aerosols. 
Also, 1.06-ym retrievals are somewhat less susceptible to transmission 
errors in the stronger layers. The principal virtue of the 0.53-ym 
wavelength is the stronger signal available (because of both atmospheric 
backscatter and quantum efficiency). These stronger signals are essen- 
tial for observing the mesospheric aerosol layers at useful horizontal 
and vertical resolutions. 
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At both 0,53 and 1,06 ym, errors in the gas density profile are the 
dominant source of uncertainty in measuring nonvolcanic stratospheric 
and upper tropospheric aerosols. Although useful retrievals of these 
aerosols can be made at 1.06 ym by using conventional density informa- 
tion, retrievals at both wavelengths would be improved by better density 
information. Part Two^ describes a method for obtaining such informa- 
tion from the lidar itself. 
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Appendix A 
ATMOSPHERIC MODELS 


A, Molecular Density 

At heights above about 30 km, our model molecular densities are 
equal to the corresponding seasonal and latitudinal models in the IJ.jS. 

on 

Standard Atmosphere Supplements , 1966 . u Below about 30 km, our model 
molecular densities are computed from temperature and pressure mea- 
surements made on each balloon flight described in Section D of this 
appendix. These balloon— measured densities differ by less than 10 per- 
cent from the corresponding seasonal and latitudinal models of the IJ.^. 

30 ~~ 

Standard Atmosphere Supplements , 1966 . 

An exception to this rule is the high-latitude case, in which the 
30- to 90-km model is for July, 60°N, whereas the balloon flight was 
made in November, at 72°N. In this case we replaced the balloon densi- 
ties with the 60°N July model densities at all heights to avoid a 
mismatch in the density profile. 


B. Ozone 


Our model ozone concentration profiles are taken from the seasonal 
and latitudinal plots assembled in the review by Wu.^ 


C. Noctilucent Cloud and Mesospheric Aerosol 

The mesospheric aerosol model is based on lidar observations made 

28 29 

at the ruby-laser wavelength 0.694 ym by Fiocco and Grams * who 
observed vertical distributions of mesospheric aerosols at high lati- 
tudes when noctilucent cloud (NLC) displays were expected. These mea- 
surements indicated that, at high latitudes, the altitude region from 60 
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to 70 km contains significant particulate material during the summer, 
and that the amount of this material appears to be slightly enhanced (by 
a factor of ~2) during periods of NLC activity. The experiments also 
provided information on the optical thickness and vertical structure of 
an NLC, 

oo 

The first series of experiment s^ 0 were carried out with a lidar 
system installed near Fairbanks and a mobile lidar system operated at 
several locations in Sweden during the summer of 1964, These experi- 
ments gave evidence for a scattering layer between 68 and 72 km on two 
consecutive nights for which NLC were observed on the horizon at each 
location. The average particulate optical thickness of the layer during 
that time was estimated to be (1 ± 0.2) x 10“ 6 . During the summer of 
1964, NLC activity was moderate, and displays were never observed over- 
head when the lidar was in operation. 

O Q 

The second series of lidar NLC experiments was conducted near 
Oslo from early July to mid-August 1966. These experiments provided 
further evidence for a broad region of enhanced backscattering 
throughout the summer in the 60- to 70-km interval with a well-defined 
maximum in the region from 66 to 68 km. The echoes from that layer were 
strongest on nights during which NLC were present — as would be expected 
from the 1964 results. The average particulate optical thickness for 
the 2-km interval from 66 to 68 km was (7 ± 1) x 10”^ for all nights 
when NLC were visually observed; this is more scattering by a factor of 
about 3 than would be expected from that interval for molecular density 
profiles specified by the Standard Atmospheric Supplements , 1966^ 

in the July standard atmosphere for 60°N. For nights without NLC 
displays, the average particulate optical thickness for that same alti- 
tude interval was estimated to be about 3 x 10“** — about one-half the 
value observed during NLC activity. The above results document the 

existence below the NLC region of a layer of particles that is normally 

32 

located near the mesopause at about 80 km. Fiocco and Grams discuss 

possible mechanisms for formation of a 70-km layer, while Grams and 
33 

Fiocco discuss the conditions necessary for formation of NLC layers at 
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the mesopause. The experimental results described above are based on an 
analysis of data obtained by summing the number of photoelectrons 
received in 2-km altitude intervals between 40 and 100 km. For the sin- 
gle case during the summer of 1966 when the lidar was operating with an 
overhead NLC display, oscilloscope photographs for each individual laser 
pulse were analyzed to provide data for a detailed study of the tran- 
sient features of the NLC, Range resolution was about 100 m. Those 
results indicated that an NLC detected by the lidar at about 75 km was 
no more than 0.5 km thick. The optical thickness of the layer was 
estimated to be 10 ^ during the short time the NLC was directly over the 
lidar system. 

We used the data in Figure 1 of Ref. 29 to derive scattering ratio 
profiles referenced to the atmospheric density profile specified in the 
60°N July model of the IJ.S^ Standard Atmosphere Supplements , 1966 . ^ Our 
"typical" high-latitude mesospheric aerosol profile, for nights without 
NLC displays, is shown in Figure A-l . Our model for nights with NLC is 
shown in Figure A-2 . In both figures, the "backscattering coefficient" 
is the radar cross section, per unit volume and 4n sr, of the aerosol 
(solid line) and the air molecules (dashed line) at the wavelength 
0 . 6 94 3 uni . 

Our choice of a size distribution for the 60- to 70-km layer was 
also guided by results presented in the same publication. Shettle and 
Fernr 5 proposed a very broad log-normal distribution function for extra- 
terrestrial meteoric debris with geometric mean radius r = 0.03 ym and 

O 

geometric standard deviation of <r = 3. If we assume the particle den- 

o 8 

sity is -2 g cm J , the results of Ref. 29 suggest that the meriodional 
circulation pattern has three effects: to remove all particles with 
radius less than ~0.02 ym; to enhance significantly the concentration of 
particles with radius between 0.02 and 0.1 ym; and to cause only minor 
changes at larger particle radii. We therefore chose a narrower distri- 
bution function with r = 0.04 ym and cr = 1.5 for both the "typical" 

o o 

and the "high-concentration" cases. For the refractive index, we used 
results for the "synthetic aerosol" model published by Ivlev and 
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90 


MESOSPHERIC REROSOL MODEL 



BACKSCATTERING COEFFICIENT (PER KM) 


FIGURE A-1 BACKSCATTERING COEFFICIENTS (PER 47T sr) AT THE WAVELENGTH 
X = 0.694 jLtm FOR THE "TYPICAL" MESOSPHERIC AEROSOL MODEL 
(CURVE WITH SOLID LINE ON THE LEFT SIDE) ALONG WITH 
COEFFICIENTS FOR MOLECULAR SCATTERING (DASHED LINE) AND 
TOTAL SCATTERING (SOLID LINE) 
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KILOMETERS 


MESOSPHERIC REROSOL MODEL 



BACKSCRTTER ING COEFFICIENT (PER KM) 


FIGURE A- 2 SAME AS FIGURE A-1 EXCEPT FOR THE NLC OR HIGH CONCENTRATION 
MESOSPHERIC AEROSOL MODEL 


3 
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3S -3b 

Popova, which was also used by Fiocco et al. to specify optical pro- 
perties of high altitude aerosols in a study of the energetic equili- 
brium of small particles in the upper atmosphere. 


For the NLC case, we added a 0.5 -km- thick layer of particles cen- 
tered at 80-km altitude to the high-concentration 70-km layer. The 
lidar backseat tering cross section of the NLC was taken to be 1.6 
x 10 5 km 1 sr 1 at \ = 0.6943 ym in accordance with the lidar NLC 
observations of Ref. 29. For the refractive index of NLC particles, we 
chose to use a real refractive index of m = 1.33 for our calculations — 
as did Grams and Fiocco JO for visible wavelength calculations in their 
study of the equilibrium temperatures of NLC particles. The NLC parti- 
cles are assumed to be somewhat larger than the extraterrestrial parti- 
cles in the 60- to 70-km layer with a larger, more sharply peaked size 
distribution, which is consistent with the expectations for particle 
formation. We used r = 0.1 pm and c r = 1.3 for the NLC particles; the 

6 o 

peak in the particle area distribution function for those log-normal 
parameters is very close to the of ten-0. 13 ym value quoted for NLC 

par tides. 

We also used our high-latitude mesospheric aerosol model to specify 
optical parameters for a midlatitude mesospheric layer. The intent here 
is to simulate particle layers near the 50-km level that have been occa- 
sionally observed by investigators using techniques other than lidar. 

We believe, in accordance with the mechanisms discussed in Ref. 32, 
that the midlatitude, 50-km layer may be a southward extension of the 
same scattering layer that Fiocco and Grams observed in the 60- to 70-km 
region at higher latitudes. In this case we would expect the layer to 
be located at a somewhat lower altitude (e.g., 50 km) at midlatitude 
locations and, furthermore, we would expect that the mixing ratio of the 
particles in the layer would be approximately the same as had been 
observed for the high- latitude layer. We therefore specify our midlati- 
tude mesospheric aerosol model by assuming that our high-latitude mesos- 
pheric layer was transported downward from the 60- to 70-km level and 
centered at the 50-km level. We furthermore assumed that this layer has 
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the same mixing ratio that it had in the higher altitude and latitude, 
i.e., that the aerosol concentration is increased in the same proportion 
as the increase in molecular density between the two altitudes. The 
proposed mechanism applies only to summertime conditions, and the molec- 
ular density profile for our midlatitude mesospheric aerosol model is 
taken from the 45°N July atmospheric model described in Ref. 30. As 
noctilucent clouds are not observed in midlatitudes, the midlatitude 
mesospheric model does not Include an NLC layer for the high- 
concentration case. 

The optical parameters used in the mesospheric aerosol model are 
summarized in Table A-l. 


D. Stratospheric and Upper Tropospheric Aerosol 


Each upper tropospheric/stratospheric aerosol model was derived 

from a two-channel vertical profile of particle number measured by a 

dustsonde (optical particle counter) flown by the University of Wyom- 
41-43 

ing. The locations and dates of the dustsonde measurements are: 


• Nonvolcanic low-latitude model: Panama (9°N), November 1973 

• Nonvolcanic mid-latitude model: Wyoming (41°N), September 1978 

• Nonvolcanic high-latitude model: Greenland (64°N), November 1978 

• Volcanic mid-latitude model: Missouri (41°N), July 1975. 


The dustsonde data are shown in Figure A-3. Each dustsonde flight 

yields vertical profiles of N 0 .i 5 and Nq. 25» where N r is the number of 

particles with radius greater than r microns. These particle number 

profiles were converted to backseat ter ing profiles at five wavelengths 

by using the optical modeling procedure described in detail by Russell 
44 

et al. This procedure divides the stratosphere and upper troposphere 
into layers, each having a range of refractive indices and relative size 
distributions. Table A-2 summarizes the three main layers (or subsets); 
between each layer is a transition layer in which optical properties are 
interpolated to yield a smooth transition. All allowed size distribu- 
tion types and refractive indices are based on previously published 
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Table A-l 


OPTICAL PARAMETERS FOR MESOSPHERIC AEROSOL MODEL: 

LOG-NORMAL DISTRIBUTION FUNCTIONS dN/d log r WITH INDICATED VALUES OF 
GEOMETRIC MEAN RADIUS r AND GEOMETRIC STANDARD DEVIATION c r 

g g 

AND INDICATED VALUES OF COMPLEX REFRACTIVE INDEX 


m 

All profiles use N^ 
observed B 


( n RE ,n IM^ * 

values that are consistent with 

values at 0.6943 urn 
n ^ 


Model 

r 8 

(ym) 


(cm 3 ) 

n RE 

n IM 

A 

(um) 

B P 

(10 i2 /km/sr) 

_ E P 

(10 ^ /km) 

B n /E n 
(sr h 

NLC 

0.1 

1.3 

1.0 

1.33 

0 

0.266 

5.16 

5.22 

0.0099 







0.355 

3.06 

2.85 

0.0107 







0.532 

2.29 

0.999 

0.0229 







0.694 

2.20 

0.458 

0.0480 







1.064 

0.912 

0.107 

0.0849 

Mesospheric 

0.04 

1.5 

1.0 

1.65 

-0.005 

0.266 

1 .48 

0.862 

0.0172 

layer 






0.355 

1.12 

0.411 

0.0273 







0.532 

0.637 

0.109 

0.0583 







0.694 

0.323 

0.0426 

0.0759 







1.064 

0.0727 

0.00934 

0.0779 
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FIGURE A- 3 DUSTSONDE DATA USED TO DERIVE STRATOSPHERIC AEROSOL 
BACKSCATTER MODELS 


Dashed lines are extrapolations. Arrows mark tropopause height. See text 
for measurement sites and dates. 









Table A- 2 


STRATOSPHERIC AND UPPER TROPOSPHERIC AEROSOL 
MODEL SIZE DISTRIBUTIONS AND REFRACTIVE INDICES 


Subset Name 

Height Range* 

Relative Size Distribution Types 

Refractive Indices 

Inner 

T + 4 km 

1 log-normal, 1 exponential, 1 truncated 

1.40-Oi, 1.42-0i 

stratospheric 

to 30 km 

power law, 3 zold, 3 modified gamma 
(each adjusted to match dustsonde 
measurement) 

1.43-0i, 1.52-0i 

Tropopause 

0.7T to 

3 segmented power law fits to 10-16 km 

1.40-0i, 1.42-0i 

layer 

T + 2 km 

data of Bigg 45,46 [41 flights] and 
Gras and Michael 4 ^ [3 flights] 

1.43-0i, 1.52-0i 

Upper 

troposphere 

0.3T to 
0.6T 

Segmented power law (Toon and Pollack 4 **) 

1.33-0i, 1.52-0i, 
1. 525-0. 005i 


T Is tropopause height. For details see Russell et al. 44 




measurements and inferences. Detailed references to the original publi- 
cations, as well as the rationale for layer height boundaries, are given 
in Ref. 44. 

As noted in Table A-2, within the "inner stratospheric" layer, size 
distributions are adjusted at each height to agree with the dustsonde 
data (Nq # i 5 and Nq. 25^ from the measurement being analyzed. This pro- 
cedure makes use of the limited information on particle size distribu- 
tion available from the dustsonde measurement, while estimating 
unmeasureed size distribution details from the models. Because the 
dustsonde measurements do not totally constrain the models, a range of 
models is allowed for each measurement height. Mie-scattering calcula- 
tions were made to transform this allowed range of models to a 
corresponding range (or set) of ratios relating backscatter to particle 
number (specifically, Nq^j^). The mean of this set of backscatter-to- 
number ratios is taken as the most probable ratio, and the standard 
deviation is taken as its lo- uncertainty. (The uncertainties were not 
explicitly used in this project.) 

The above process was performed on the dustsonde data of Figure A-3 
for each of the wavelengths (1.064, 0.694, 0.532, 0.355, 0.266 ym) used 
in this study. Figure A-4 shows the resulting multiwavelength back- 
scatter profiles. The upper stratosphere — above 40 km — is not included 
because of the lack of aerosol measurements there. Between 40 km and 
the top of the dustsonde data (22 to 30 km, depending on the particular 
measurement), the data have been extrapolated. 

As a test of the optical modeling procedure, we have compared the 
derived backscatter profiles at wavelength 0.694 ym with ruby lidar mea- 
surements made at the time and place of the dustsonde measurements 
(except for the Panama measurement, for which no lidar measurements were 
made). In each case, the measured and modeled scattering ratio profiles 
agreed to within the combined measurement and conversion uncertainties 
(roughly speaking, an uncertainty of 0.03 in scattering ratio). 
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BACKSCATTER MIXING RATIO 


FIGURE A-4 MODEL BACKSCATTER MIXING RATIO PROFILES DERIVED FROM 
THE DUSTSONDE DATA OF FIGURE A-3 AND THE OPTICAL MODEL 
OF TABLE A- 2 
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Although this comparison does not prove the validity of the backscatter- 
ing profiles derived at other wavelengths, it lends confidence in the 
modeling procedure and is the best test than can be made on the basis of 
currently available data. 


E. Saharan Aerosol 

The Saharan aerosol layer consists of crustal aerosol particles 
that are generated in the deserts of northwest Africa and transported 
across the tropical Atlantic ocean north of the equator. The Saharan 
aerosol layer is generally confined between the tradewind inversion at 
~1 km and an altitude of 4 to 6 km, with a well-defined top. The tran- 
sport is most important during summer. 

The Saharan aerosol model is based on a series of measurements of 
the refractive indices, size distribution, concentrations, and vertical 
profiles of the Saharan aerosol layer that have been made in the eastern 
Atlantic (Sal Island and the GATE study area) and in the western Atlan- 
tic (Barbados and the BOMEX study area). From these data, we consider 
three cases of aerosol concentration: high, typical, and low. The 

molecular density profiles are assumed to be given by the 15° annual 

30 

model of the atmosphere. 


The size distributions in each case have been chosen to be single 
modal log-normal distributions of the form 


dN 


d log r 


N t 


yflvt log (Tp 


exp 


(log r - log r g ) 


2 log <r. 


8 


(A-l) 


with N, rg, and o*g being, respectively, the number (per cm ), the 
geometric mean radius (^m) , and the geometric standard deviation. The 
parameters of the distributions are shown in Table A-3. The size dis- 
tribution for the typical concentration was chosen on the basis of vari- 

49 

ous continental measurements reported by Patterson and Gillette as 
well as the series of measurements during GATE reported by Savoie and 
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Table A- 3 


SAHARAN AEROSOL MODEL 



r g 

(p») 

°k 

V 3 

(cm 3 ) 

X 

(pm) 

n R 

n IM 

Low concentration 

0.30 

2.00 

1.00 

1.064 

1.530 






0.694 

1.550 






0.532 

1.560 






0.500 

1.560 






0.355 

1.580 

-0.022 





0.266 


-0.048 

Typical concentration 

0.35 

2.20 

1.00 

1.064 

1.530 






0.694 

1.550 






0.532 

1.560 

-0.006 





0.500 

1.560 






0.355 

1.580 

-0.022 





0.266 


-0.048 

1 

High concentration 

0.60 

2.10 

1.00 

1.064 

1.530 

-0.004 





0.694 

1.550 

B i 

i 




0.532 

1.560 






0.500 

1.560 

Cfjl 





0.355 

1.580 

-0.022 





0.266 

1.600 

-0.048 


Nominal value of particle concentration; actual values are 


0un -1 8r -1 ) 


-1 

(k® *> 


0.201E-02 

0.195E-02 

0.185E-02 

0.181E-02 

0.168E-02 

0.161E-02 

0.347E-02 

0.325E-02 

0.308E-02 

0.304E-02 

0.289E-02 

0.283E-02 

0.831E-02 

0.762E-02 

0.734E-02 

0.731E-02 

0.717E-02 


(km l ) 

0. 193E-02 
0.183E-02 
0.164E-02 
0.156E-02 
0.113E-02 
0.888E-03 

0.326E-02 
0.296E-02 
0.261E-02 
0.245E-02 
0. 179E-02 
0.154E-02 

0.760E-02 

0.667E-02 

0.587E-02 

0.549E-02 

0.422E-02 


0.1 14E-03 
0.17 IE-03 
0. 169E-03 
0. 148E-03 
0.651E-04 
0. 152E-04 

0.235E-03 
0.278E-03 
0.237E-03 
0, 198E-03 
0.705E-04 
0. 147E-04 

0.610E-03 

0.596E-03 

0.429E-03 

0.337E-03 

0.790E-04 


0.567E-01 
0 .879E-01 
0.916E-01 
0.816E-01 
0.388E-01 
0 • 942E-02 

0.678E-01 

0.858E-01 

0.771E-01 

0.650E-01 

0.244E-01 

0.521E-02 

0.733E-01 
0.783E-01 
0.585E-01 
0.461E-01 
0.1 10E-01 


0.702E-02 


0.380E-02 


0. 167E-04 


0.238E-02 
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Prospero. The distribution for our typical concentration is approxi- 
mately that reported by Savoie and Prospero as the geometric mean for 

the distributions observed at Sal Island, These data, as well as the 

51 

results discussed by Grams et al., indicate that the highest concen- 
trations are characterized by the presence of relatively large parti- 
cles, so our high concentration model has a slightly larger mass mean 
radius. The size distribution for the low concentration model is 
representative of the conditions measured on Barbados by Savoie and 
Prospero; the characteristic sizes are smaller than the other cases, 
reflecting losses of larger particles during their transport across the 
Atlantic . 


The refractive indices of the Saharan aerosol model are those 

52 53 49 

reported by Patterson et al . and by Patterson. Patterson et al, 

reported the real index of refraction n R g at two wavelengths: 0.63 pm 
and a rather broad interval centered at 0.55 pm. From these data 
points, a wavelength dependence was estimated for the interval between 
0.45 and 0.7 pm. This estimated n RE was used for the laser wavelengths 
of 0.53 and 0.69 pm; the modeled n RE at 0.266 and 0.355; and 1.06 pm was 
determined by extrapolation. The imaginary index of refraction n IM was 
measured between 0.25 and 0.70 pm as shown in Figure A-5. nj^ at 1.06 
pm was assumed to be equal to that at 0.7 pm, an assumption based on 
several measurements that show little variation in njji between 0.7 and 
1.0 pm (Ref. 54, for example). Measurements of Saharan aerosol samples 
collected on both sides of the Atlantic showed no significant variation 
in either real or imaginary index with distance from the source; conse- 
quently, we chose the same wavelength-dependent refractive indices for 
our models. Also, since the crustal aerosol is nonhygroscopic , no 
corrections to these measured optical properties to account for varia- 
tions in relative humidity were applied. 

Although the size and refractive index measurements discussed above 
show a great deal of uniformity, the concentrations and vertical dis- 
trubtions of the Saharan aerosol can vary widely. Some measured rela- 
tive concentration profiles for the eastern Atlantic and the western 
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IMAGINARY INDEX OF REFRACTION 
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FIGURE A- 5 n |M VALUES DETERMINED FOR CRUSTAL AEROSOLS 
FROM VARIOUS GLOBAL SOURCES 

Saharan aerosols ( , 0); North American aerosols ( , x); 

and a southwest Asian desert aerosol ( ). Also shown is a relative 

n, M curve determined from absorption data (■••) (from Patterson 53 ). 


800 
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SAHARAN AEROSOL- EASTERN ATLANTIC 



RELATIVE CONCENTRATION 

(a) 


SAHARAN AEROSOL- WESTERN ATLANTIC 



RELATIVE CONCENTRATION 

(b) 


FIGURE A-6 RELATIVE CONCENTRATION PROFILES OF THE SAHARAN AEROSOL 

MEASURED IN THE EASTERN ATLANTIC (a) AND IN THE WESTERN ATLANTIC (b) 

The data are those of Grams et al. 51 (•■•, — ), Kondratyev et al. 55 ( ), 

and Prospero and Carlson 56 ( , ). 


Ob 

Ui 





Atlantic are shown in Figure A-6(a,b). The profiles in Figure A-6(a) 

are derived from Kondratyev"^ and Grams (unpublished data); those in 

56 

Figure A-6(b) were reported by Prospero and Carlson, The Kondratyev 
profile is based on a measurement of relative number, that of Grams on 
relative light scattering cross sections, and that of Prospero on rela- 
tive mass concentrations. Each profile is consistent, having a base at 
approximately 1 km and a top at 4 to 5 km. The layering in these pro- 
files is a common feature of the aerosol distribution, but the details 

are quite variable. For our model, we chose the profile published by 
55 

Kondratyev for the high and the typical concentration cases and the 

56 

layered profile published by Prospero for the low concentrations 

cases. The aerosol concentrations in our model were normalized to the 

total optical depth of the Saharan layer, as measured during summer. We 

chose optical depths of approximately 1 for high concentrations, 0.5 for 

typical concentrations, and 0.2 for low concentrations. The high con- 

57 

centration value is based on the averages of Carlson and Caverly for 
high turbidity days at Sal Island during GATE; the typical concentration 
optical depth is based on average turbidities at Sal Island during GATE; 
and the low concentration value is based on average turbidities for Bar- 
bados during GATE. These summer measurements represent maximum turbidi- 
ties for the time at which the transport of dust is at a maximum; winter 
measurements may be lower by as much as a factor of ten. In addition, a 
time series of concentration measurements for the Barbados samples for 
the years 1965-75 shows that the GATE measurements of summer 1974 were 
among the highest in the time considered. The lowest summer mea- 
surements were lower by approximately a factor of 3 than the summer max- 
ima considered in our model. The winter minima were approximately equal 
in all of the years. 


F. Marine Aerosol 

The size distributions and refractive indices used in the marine 

48 

aerosol models are based on the work of Toon and Pollack and Shettle 
34 

and Fenn. Table A-4 summarizes these properties. 
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Size Distribution 


Refractive Index 


B p ( A)/E p (X) - 


A (pm) - 


Size Distribution 

Refractive Index 

0.266 

0.353 

0.532 

Segmented power 

law [Eqs. (A2,A3)]* 

1.52-Oi 

1.525-0.0051 

0.069 

0.037 

0.067 

0.038 

0.067 

0.041 

Biaodal log 

normal [Eq. ( A4 ) ] ^ 

1.52-Oi 

1.525-0.0051 

0.062 

0.031 

0.065 

0.035 

0.072 

0.043 


Weighted average 

0.055 

0.056 

0.060 


Standard deviation, 
X 

3!. 

28. 

24. 
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The equation of the Toon and Pollack size distribution is 

C r < 0 . 1 ym 


n(r) = 


vm ) 


\ v+1 

C | ym} 50 ym - r 0.1 ym 


(A- 2) 


0 


50 ym - r 50 ym. 


where 


2.0 0.1 p < r <5,0 p 

4.0 5.0 ym r <_ 50 ym * (A-3) 


and n(r)dr is the number of particles with radii between r and r + dr. 
48 

Toon and Pollack developed this equation as a best-fit solution to a 
collection of tropospheric aerosol measurements selected to represent 
the global background aerosol for heights below 3 km. Because most of 
the Earth's surface is covered by oceans, we adopted this distribution 
as a likely descriptor of the marine aerosol (i.e., the aerosol present 
over the ocean surface, whether of continental or marine origin). 


34 

The equation of the Shettle and Fenn marine aerosol is 


n(r) = 2 


Ni 


exp - 


i=1 In 10 ro*i 


(log r - log r^) 
— 

2<r ± 


(A-4) 


where 

-4 

^2 = Nj x 4.71 x 10 
rj = 0.005 ym 
r 2 “ 0.3 ym 
= 0.475 
= 0.4. 

The small-particle mode [i = 1 in Eq. (A-4)] represents particles of 
continental origin; the large-particle mode [i = 2] represents those of 
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sea-spray origin. Although this suggests that the two modes are likely 
to have different refractive indices, for the purposes of our model we 
keep refractive index independent of particle size. 

48 

With guidance from Toon and Pollack we modeled the marine aerosol 

composition as one-third soil particles (taken as basalt), two-thirds 

sea salt (taken as NaCl), and traces of sulfate [taken as (^ 4 ) 2804 ]. 

The refractive index of the soil particles is taken as 1.525- 
48,58,59 

0.005i and that of the sea salt and sulfate as 1.52-0i (see 

Table A-4). The latter value is somewhat less than that quoted in the 
literature for pure NaCl (1.58-1.53 — see Toon et al. ) and allows for a 
small water admixture caused by the deliquescence of NaCl and the typi- 
cally high relative humidity of marine environments. We modeled these 

refractive indices independently of wavelength. The data of Toon et 
60 

al . show that, for NaCl and (NH 4 ) 2 S 04 , this approximation is good to 
an accuracy of ± 0.02 for the real part and 10 for the imaginary part 
over the wavelength region 0.35 to 1.06 ym. 

Mi e- scattering calculations were made for each combination of size 
distribution and refractive index shown in Table A-4. Each combination 
yielded the multiwavelength backscatter-to-extinction and backscatter- 
to-number ratios shown. Weighted averages of these ratios were then 
taken and used to derive multiwavelength backscattering profiles from 
specified midvisible optical thickness values and relative vertical pro- 
files, as described below. Our reason for averaging the results of the 
Mie-scattering calculations, rather than averaging the refractive 
indices input to the Mie calculations, was to avoid errors caused by the 
nonlinear dependence of scattering and extinction on refractive index. 
These errors are discussed at some length by Bergstrom and Toon et 



The vertical distribution of our marine aerosol models was based on 

62 

the lidar measurements reported by Livingston and Uthe, who noted that 
the surface-based marine haze layer in the Atlantic and Mediterranean 
had a mean height of 0.5 km with a standard deviation of 0.3 km for the 
22-day cruise they reported. They also noted that elevated layers below 
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3.6 km in altitude were usually present, and that as a result the 
surface-based layer typically accounted for less than one-half of the 
vertically integrated backscatter below 3.6 km. 

With these results in mind, we adopted a three-layer model for the 
marine aerosol. As shown in Figure 2 of the main text, the surface- 
based layer has constant backscatter mixing ratio up to a height of 0.4 
km for the high-latitude version and 0.6 km for the mid- and low- 
latitude versions. The elevated layers are independent of latitude and 
centered at heights of 1.5 and 2.5 km, with approximately Gaussian 
shapes and half-width (HWHM) of 0.5 km. These layers have equal peak 
values for backscatter mixing ratio, with the peak value adjusted so 
that the surface-based layer contains only about 50 percent of column 
backscatter below 3.6 km. 


The absolute backscattering coefficients in our marine aerosol 
models were determined by combining backscatter-to-extinction ratios for 
the averaged model optical properties (Table A-4) with the requirement 
that the midvisible (A = 0.53 pm) optical depth of the marine aerosol be 
0.16 for the high-latitude case and 0.19 for the low- and mid-latitude 
cases. These values may be compared with a large set of time-averaged 

maritime tropospheric midvisible optical depths summarized by Toon and 

48 48 

Pollack. The values selected by Toon and Pollack to represent cases 

with little continental influence ranged from 0.02 to 0.17, with the 
larger values tending to occur at low latitudes. The midvisible parti- 
culate extinction coefficient E (0.53 ym) at the surface for each of our 

y -1 

marine aerosol models is 0.167 km . Using the fact that Ep(0.53 ym) 

» E p (0.55 ym) and that the gaseous Rayleigh extinction coefficient 
E S g(0.55 ym) sj 0.012 yields a total extinction coeficient. 


E(0.55 ym) « 0.18 km 


-1 


(A-5) 


90 



V 


3.91 


The Koschmleder relation. 


E(0.55 ym) 


(A-6) 


yields a surface meteorological range V of 22 km (13.5 miles) for each 
model. As noted by Elterman and Collis and Russell, such a meteoro- 
logical range is associated with very light haze conditions, usually 
described as "clear." Thus, although our marine aerosol models 

tend to have somewhat larger optical thicknesses than the globally aver- 

48 

aged, no-continental-influence values of Toon and Pollack, they 
nevertheless represent light haze conditions. 


G. Cirrus Clouds 

In contrast to our aerosol models, our cirrus cloud models are 
based on empirical (visual, lidar, radiometric, in-si tu) observations of 
cloud macroproperties , rather than on calculations from cloud m icropro- 
perties (ice crystal number density, size distribution, shape, and 
orientation). Since clouds are highly variable, we have included three 
cirrus layers of different density, height, and geometrical thickness in 
each latitude band of the model, rather than trying to describe some 
"typical" state. Figure A-7 shows the cirrus layers in terms of mul- 
tiwavelength backscatter mixing ratios, each superimposed on the nonvol- 
canic aerosol model for the corresponding latitude band. Table A-5 
lists several properties of these layers. 

The properties listed in Table A-5 were derived as follows. Each 
layer's top height is given by h^, and its geometrical thickness by 
The mean midvisible gas backscattering coefficient in the layer, B g (0.53 
ym), is computed from the corresponding profile of gas density D using 

Bg(0.53 ym) = y( 0.53 ym) P(h^»>Ah) , (A-7) 
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ALTITUDE 


vO 

N3 



BACKSCATTER MIXING RATIO. 


FIGURE A- 7 CIRRUS CLOUD MODEL BACKSCATTER MIXING RATIOS SUPERIMPOSED 
ON THE MODEL AEROSOL IN EACH LATITUDE BAND 






Table A-5 


PROPERTIES OF MODEL CIRRUS CLOUD LAYERS 


hfp 

(km) 

Ah 

(km) 

R(0.53 pm) 

B (0.53 pm) 
Sr 1 sr' 1 ) 

(m 1 sr 

r* 

r+ 

Viability* 

w s 

(g rn 3 ) 

High-Latitude 

9 

1 

1.16 

6 x 10“ 7 

1 x 10 -7 

0.005 

0.002 

s 

1.5 x 10 -4 

7 

2 

1.67 

8 x 10 -7 

5 x 10 -7 

0.05 

0.02 

s 

7 x 10 -4 

4 

1 

11.00 

1 x 10~ 6 

1 x 10“ 5 

0.5 

0.2 

v 

1.5 x 10 -1 

Mid-Latitude 

12 

1 

1.16 

4 x 10~ 7 

5 x 10“ 8 

0.013 

0.0012 

s 

7 x 10 -5 

10 

2 

1.67 

6 x 10 -7 

4 x 10 -7 

0.04 

0.016 

s 

6 x 10 -4 

7 

3 

11.00 

8 x 10“ 7 

8 x 15 

1.2 

0.48 

V 

1 x 10" 1 

Low-Latitude 

16 

1 

1.16 

1 x 10 -7 

2 x 10 -8 

0.001 

0.0004 

s 

5 x 10' 5 

14 

2 

1.67 

3 x 10' 7 

2 x 10“ 7 

0.02 

0.008 

s 

3 x 10" 4 

11 

3 

11.00 

6 x 10 -7 

6 x 10" 6 

0.9 

0.36 

V 

2 x 10 -1 


Assumes intrinsic backscatter-to-extinction ratio, T “ 0.02 sr -1 . 

^ Assumes effective backscatter-to-extinction ratio, T’ “ 0.05 sr 
z 

Based on Bp, ^h, and previous comparisons of lidar, human observer, 
and satellite TV data. See Table A-6. 

8 

Assumes effective particle radius, a ■ 50 pm. 


where Y is the Rayleigh backseat ter-to-mass ratio for air, given by: 

Y(0 . 53 pm) = 1.30 x 10 9 m 2 sr 1 g 1 . (A-8) 

Cloud particle backscattering coefficients are then derived using the 
specified scattering ratios R and 

Bp = Bg(0.53 pm) x [R(0.53 pm) - 1] . (A-9) 
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Note that we have omitted a wavelength descriptor from the left side of 
Eq. (A-9) because the size of cirrus particles is large compared to the 
wavelengths used in this study. We therefore used the geometric optics 
approximation that cloud particle backscatter and extinction were 
independent of wavelength. 

Optical thicknesses t are obtained from 

T = Bp Ah/r , (A-10) 

where P is the cloud particle backscatter-to-extinction ratio. Follow- 

65 

ing the empirical data of Platt, we modeled cirrus clouds as having 

P = 0.020 sr" 1 . (A-ll) 


Again, because cirrus crystals are large in comparison to the 

wavelengths used in this study, lidar measurements of cirrus can be 

strongly affected by multiple scattering. That is, light-scattering by 

cirrus particles is very strongly forward-peaked, thus returning much 

scattered light to (or near) the original direction of travel. These 

multiple scattering effects can be described in an approximate way by 

using an effective backscatter-to-ext inction ratio P' in computing lidar 
" ^ 

pulse transmissivities. Platt has had success using 

P' =0.05 sr” 1 ( A-l 2 ) 

in his studies, and we have adopted his value in our models. Note that 
P' exceeds P because multiple scattering tends to decrease the effective 
extinction coefficient for a given backscatter coefficient. Lidar pulse 
transmission is described by an effective optical thickness, defined 
analogously to Eq. (A-10) as: 

x'=B p Ah /r . (A-l 3) 
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To aid in relating our model cloud layers to observable physical 
characteristics, the last two columns of Table A-5 give approximate 
descriptors of cloud visibility and water content. A cloud's visibility 
to ground-based human observers can be related to backscatter or optical 
thickness values only very approximately. Nevertheless, some useful 
distinctions can be developed on the basis of the large set of lidar, 
human observer, and satellite TV observations acquired in many seasons 
and latitudes by Evans. These distinctions are summarized in Table A-6. 
On the basis of Table A-6 we classified each of our model cirrus layers 
as either visible or subvisible (to a ground-based human observer). As 
can be seen from Table A-5, in each latitude band of the model the top 
two cloud layers are nominally subvisible, whereas the bottom layer is 
nominally visible. Further, the top layer in each band is similar in 
optical and physical density to the weakest cloud described by Uthe and 
Russell. However, the model layer is geometrically thicker than the 
measured layer. 


Cloud water content can be approximately derived from the back- 
scatter values by treating the particles as spheres and assuming an 
effective particle radius, a. Again, using the geometric optics approx- 
imation, one can relate particle number density N and extinction coeffi- 
cient Ep by 



2na 


(A-14) 


Cloud water content is given by 


W = — TTa pN 
3 


(A-15) 
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Table A-6 


APPROXIMATE RELATIONSHIPS BETWEEN CIRRUS CLOUD BACKSCATTER, 
THICKNESS, AND VISIBILITY* 


Particle Backscatter Coefficient (m * sr *) 

i— > 

o 

i 

■^j 

io " 6 

10 5 

10 -4 

-3 

10 


Weak 

Cirrus 

Medium 

Cirrus 

Strong 

Cirrus 


Never 

observed 

visually 

Daytime 
visible 
only if 
Ah > 2 km 

Daytime 

visible 

if 

Ah > 0.5 km 

Always daytime 
visible; may even 
appear opaque 



Seldom Moonlight visible 

detected under good 

via conditions if 

satellite TV Ah > 1 km 


Based on Reference 1. 


where p is particle mass density (0.92 g cm for ice, 1 g cm for 
water). Combining Eqs. (A-14) and (A-15) with the definition of the 
backscatter-to-extinction ratio. 


r s b p /e p , 


then yields 


2 ap 

W B p . 

3 r 


(A-16) 


(A-17) 


For ice crystals, this in turn yields 

-3 -1 -1 

W(g m ) = 30a(ym)Bp(m sr ) . (A-18) 

The W values shown in Table A-5 were derived from Eq. (A-18) using an 
effective particle radius a of 50 um. By comparison, the in-situ 
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r 


measurements of Varley and Brooks showed that the particle radius con- 
tributing most effectively to ice content in sampled cirrus clouds 
ranged approximately from 10 to 100 ym. 
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Appendix B 

BACKS CATTER-TO-EXTINCTION RATIOS 


The lidar data analysis method assumed in this study constructs 
particulate extinction profiles E p (z) by dividing lidar-derived back- 
scatter profiles B p (z) by a possibly height-dependent model 
backscatter-to-extinction ratio, Y p (z). Thus, the uncertainty in Y p 
contributes to the uncertainty in particulate optical depth x • In 
fact, for the cases of practical concern in this study, where 0.01 _< x p 
1.0, S^p/Yp is the dominant contributor to 6x p /xp. In some cases this 
dominance results from a need to solve iteratively for B p (z) and E p (z) 
profiles that are consistent with the model Y p ( z ) profile. 

Because of the importance of 5Y p /Y p in determining 6x p /T p [and 
hence the transmission-induced uncertainty in B p (z) — see Eqs. (12), 

(20), and (24)] we have surveyed model and measured values of Y p a s a 
means of evaluating the possible range of values (and hence the uncer- 
tainty in any chosen model value). Also, we have searched for correla- 
tions between Y p and the parameter a defined by 

-a 

B p (0. 53 ym) /o.53\ 

B p ( 1.06 ym) \1.06/ (B— 1) 

or 

ln[Bp(0.53 ym)/Bp(1.06 yin)] 

a = - . 

In 2 ( B " 2 > 

Our reason for investigating this correlation is that a can be deter- 
mined from a lidar measurement of B p (0.53 ym)/B p (1.06 ym); hence, if a 
and Yp are correlated, knowledge of B p (0.53 ym)/B p (1.06 ym) can be used 
to reduce uncertain! ty in Y p . [For an illustration of this process, see 
Ref. 44.] 
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The results of this survey are shown in Figure B-l • Figures B-l(a) 
and B-l (c) show model values of Y p (0.53 pm) and Y p (1.06 ym) calculated 
for the various aerosol and cloud models used in this study, each plot- 
ted versus the corresponding value of a . (For details on model size 
distributions and refractive indices, see Appendix A, Sections A-3 
through A-7.) These results indicate that Y p does not depend very 
strongly on wavelength. Figure B-l (b) shows values of ^(0.69 ym) mea- 
sured in experiments using ruby lidars. Since a was not measured in 
these experiments, Y p (0.69 ym) cannot be plotted versus a, but neverthe- 
less the range of Y p (0.69 ym) can be seen. 

A. Stratospheric Aerosols 

Many comparison studies have shown that stratospheric aerosols are 
described fairly well by spherical-particle models of the type used to 
construct the curves in Figure B-l(a,c). (For a summary of such com- 
parisons, see Appendix A and Ref. 44.) Thus, although measurements of 
stratospheric Y p are not available, the model results in Figures B- 
l(a,c) can be adopted with some confidence. These results indicate that 
for the stratospheric aerosol, cl is a good predictor of Y p (0.53 ym) , 
provided a 1.5, as it is for most stratospheric aerosols. However, 
the practical value of this predictive ability is diminished by the 
facts that: 

• For the background stratospheric aerosol, a is difficult to mea- 
sure by lidar because B p (0.53 ym) is difficult to measure. 
[Uncertainties in molecular density severely degrade mea- 
surements of B p (0.53ym) at most stratospheric heights, even when 
using the three-wavelength analysis technique.] 

• Highly accurate particulate backscatter-to-extinction ratios are 
not required for deriving lidar backscatter profiles in the 
unperturbed stratosphere, because stratospheric particulate 
transmission is nearly unity and thus is a small source of 
uncertainty in analyzing lidar data. 
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FIGURE B-1 MODEL AND MEASURED VALUES OF PARTICULATE 
BACKSCATTER-TO-EXTI NOTION RATIO, 7 

'p 

(a) Model values for X - 0.53 (xm, plotted versus a defined by Eq. (B-2); (b) Measured 
values for \ = 0.69 jum. Letters indicate locations (all lower tropospheric). A = Tucson, 
Arizona (Refs. 68-71); S = San Francisco (Ref. 20); W = Seattle, Washington (Ref. 73); 

E = Hempstead, England (Ref. 74); (c) as in (a) but for A = 1.06 fxm. 


B. Tropospheric Aerosols 

The most critical need for a well-defined relationship of Yp versus 
a is in the troposphere, because it is here that aerosol transmission 
uncertainties become significant error sources in analyzing lidar data. 
Here, a can frequently be measured more accurately than in the strato- 
sphere. Unfortunately, Figures B-l(a,c) show that, for the tropospheric 
aerosol models, the correlation between Yp and a is very weak. In fact, 
any systematic variation of Yp appears to be as well predicted by loca- 
tion (plus backscatter strength) as by a — that is, the marine and 
Saharan values for Yp systematically exceed those for the upper tropo- 
sphere (excluding the Saharan cases). The reason that there is no 
well-defined relationship for Yp versus a within any of these groups is 
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the large range and sometimes strong wavelength dependence of refractive 
index (real and imaginary) in the troposphere as compared to the strato- 
sphere . 

Given this result, one can adopt the alternate approach of select- 
ing Yp on the basis of aerosol location, independent of a . A measure of 
the uncertainty for this approach can be obtained by calculating the 
mean and standard deviation of the appropriate groups of data points in 
Figures B-l(a,c). The results, shown in Table B-l(a,b) yield lo* uncer- 
tainties ranging from 19 to 37 percent of the mean. 

However, this approach may underestimate tropospheric 6T p/7p, and 
overestimate Y p , because the input set of models is limited — in particu- 
lar by excluding all highly absorbing urban aerosols and all nonspheri- 
cal aerosol particles. Such particles typically have smaller Y p values 
than do spherical-particle models with little or moderate absorption. 
This concern is supported by the measured lower tropospheric data points 
in Figure B-l(b), which are systematically smaller than the correspond- 
ing model points in Figures B-l(a,c). It would be useful to calculate 
the mean and standard deviation of all the measured data points in the 
experiments referred to in Figure B-l(b). However, this is not possible 
because many of the references do not give individual data points. An 
approximate approach is to calculate the mean and standard deviation of 
two sets shown in Figure B-l(b): (1) The eight "central" data points; 
and (2) the central data points plus the ten upper and lower error bar 
extremes. The results are shown in Table B-l(c). Note the smaller 
values for Y p and larger values for 6Y p /Y p , as compared to the model 
values in Table B-l(b). 

C. Clouds 


The only cloud value for Y p included in Figures B-l(a,c) is the 

cirrus value 0.05 sr ^ (independent of wavelength). This is an effec- 

65 

tive value adopted by Platt and by the authors of this study to 

74 75 

include multiple scattering effects. Recently, several studies * 
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Table B-l 


MODEL AND MEASURED AEROSOL PARTICLE BACKSCATTER-TO-EXTINCTION RATIOS 


(a) 

Model Values, Upper Troposphere (excluding Saharan) 
[Figs. B1 (a ,c) ] 

X(ym) 

ysr" 1 ) 

%<sr -1 ) 

4 yy*> 

0.53 

0.023 

0.007 

31 

1.06 

0.023 

0.005 

19 


(b) Model 

Values, Marine and Saharan 



[Figs. Bl(a ,c) ] 


X(ym) 

ysr- 1 ) 

^YpCsr -1 ) 

6Y p /y%) 

0.53 

0.054 

0.020 

37 

1.06 

0.053 

0.016 

32 


(c) Measured Values, Lower Troposphere, 

— 0.69 ym 



[Fig. Bl(b) ] 






_ * 

Set 

Y ( s r l ) 

<5Y (sr A ) 

5Y / Y (%) 


P 

P 

P P 

i 

0.028 

0.014 

50 

2 

0.033 

0.021 

66 


See text for definition of sets. 


have shown that single-scattering Y p values for a wide range of spheri- 
cal water-droplet cloud and fog models vary by only about 15 percent 
around the mean value 0.06 sr"*^ for wavelengths between 0.5 and 1.06 ym. 
However, this small variation must be viewed with caution in lidar stu- 
dies. First, water-droplet clouds are usually optically dense, so that 
multiple scattering effects must be taken into account in attempting to 
infer quantitative values of B p from the lidar return. Second, although 
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cirrus clouds are frequently tenuous, the large size of the crystals 
produces very strong forward scattering (a multiple scattering effect 
that decreases effective extinction); also, the crystals can be markedly 
aspherical. Thus, the small variation of single-scattering spherical 
water-droplet Yp is not applicable to the lidar measurements in this 
study. Given the wide range of crystal shapes, orientations, and sizes, 
as well as cloud optical densities, it appears that the relative uncer- 
tainty in cirrus Yp must be considerably larger than 15 percent. 


D. Summary 

On the basis of the above results, we have adapted 50 percent as a 
rough measure of the relative lcr uncertainty in Yp, and hence in Tp, for 
our simulations. (See the first paragraph of this appendix.) However, 
it can be seen that this value should be revised upward or downward in 
particular circumstances, and hence this is an area for careful atten- 
tion in individual measurements. 
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Part 3 

ORBITING LIDAR SIMULATIONS: 

DENSITY, TEMPERATURE, AEROSOL AND CLOUD MEASUREMENTS 
BY A WAVELENGTH-COMBINING TECHNIQUE 



I INTRODUCTION 


Part Two of this report simulated aerosol and cloud retrievals for 
a proposed space shuttle lidar using independent-wavelength analysis 
techniques and conventional density information. It showed that for 
upper tropospheric and nonvolcanic stratospheric aerosols in areas 
remote from radiosonde soundings (i.e., over most of the globe), errors 
in the density profiles inferred from gridded meteorological analysis 
are usually a leading source of error in retrieved particle backscatter 
profiles. Moreover, obtaining a density profile from gridded meteoro- 
logical analysis for each lidar measurement location and time entails 
considerable expense and delay. The alternative — using model density 
profiles — reduces expense and delay, but even if the models are care- 
fully selected by location and season, the resulting errors exceed those 
of gridded analysis. Clearly, the quality, speed, and cost- 
effectiveness of retrieving cloud and aerosol information from lidar 
measurements would be enhanced by density information from the lidar 
itself, provided the lidar density errors were less than those from 
models or gridded analysis. 

Previous experimental studies* ^ demonstrated that single- 
wavelength, low-spectral-resolution lidar measurements can yield very 
useful relative density profiles and temperature profiles in regions of 
the atmosphere where particulate backscattering at the lidar wavelength 
is negligible compared to gas backscattering. However, these single- 
wavelength techniques clearly fail in obtaining density or temperature 
profiles in regions where particulate backscattering at the laser 
wavelength is a nonnegligible fraction of gas backscattering, and they 
have no means of showing when this particulate contamination occurs. 

This shortcoming can jeopardize the credibility of any unusual results, 
which may be of the greatest interest (e.g., unusual density striations 
or wavelike structures). 
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DeLuisi et al • presented a technique for combining lidar signals 

at three wavelengths and two zenith (or nadir) angles, so as to obtain 

separate profiles of molecular density and particulate backscatter and 

extinction. They also showed how the molecular density profile could be 

combined with the hydrostatic equation and the ideal gas law to obtain a 

5 

temperature profile (as in Ref. 3). The DeLuisi technique assumes a 
specific one-parameter form (linear for the published cases) for the 
wavelength dependence of particulate backscattering , and uses two of the 
three lidar wavelengths to solve for the wavelength-dependence parame- 
ter. The effect of actual aerosol populations that do not follow the 
assumed functional fora for wavelength dependence was not investigated. 
Thus, the size of density and temperature errors caused by various aero- 
sol concentrations, wavelength dependences, and lidar parameters could 
not be evaluated. 

This paper presents an alternate, but related, technique for 
obtaining profiles of gas density and temperature, as well as aerosol 
and cloud backscattering. A major goal of this technique is to make the 
minimum possible extension to single-wavelength, single-angle analysis 
techniques that have been applied for many years, and for which the 
practical error sources are fairly well understood. As will be shown, 
this technique also produces large density and temperature errors when 
aerosol or cloud concentrations are sufficiently large. However, it 
minimizes these errors by using a short laser wavelength (to strengthen 
gas backscatter), and it automatically provides the user with an esti- 
mate of these errors at all heights by using cloud and aerosol informa- 
tion determined from a simultaneous longer-wavelength measurement. (It 
also reduces errors somewhat by using the longer-wavelength measurement 
to correct for short-wavelength aerosol effects.) 

After describing the technique we present an error analysis that 
includes the effects of errors in lidar signal profiles, transmission 
profiles, density normalization, and reference pressure or temperature, 
as well as aerosol or cloud contamination of the lidar density profile. 
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Finally, we simulate the measurement and retrieval process in order to 
test the algebraic error expressions and indicate the range of results 
possible. 


For purposes of illustration we assume 1.064 and 0.355 ym as the 
primary wavelengths for input to the retrieval calculations. These 
wavelengths are of special interest because: 

• They are produced by the Nd:YAG laser (plus frequency-tripling 
optics), a source that is becoming increasingly popular for 
current and proposed lidar systems (ground-, air-, and 
space-based ) J 

• They are sufficiently separated that gas backscattering is much 
stronger at one than the other. 

• The tripled wavelength, 0.355 ym, while short enough to be 
strongly backscattered by gas molecules, is just long enough to 
avoid significant absorption by the Hartley-Hug gins bands of 
ozone. 

However, the analysis techniques are quite general and could be applied 
to any other set of wavelengths with some or all of these attributes. 
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II SOLUTION TECHNIQUE 


A. Gas Density and Particle Backscatter Profiles 

Figure 1 illustrates the solution technique. The numbers shown 
indicate the following steps: 

Step 1 : Use the long-wavelength signal profile S(X^ 5 z) and a conven- 

tional density profile p c (z) (model or interpolated from mea- 
surements) to solve for the long-wavelength scattering ratio 
profile R(A^,z), and its uncertainty <5R(A,z), using the con- 
ventional single-wavelength technique [Eqs. (4) to (9) and 
(11) of Part Two of this report]. 

Step 2 } Use the long-wavelength scattering ratio profile and uncer- 
tainty with an assumed, typical wavelength dependence (and 
uncertainty) of particulate backscattering to estimate the 
short-wavelength scattering ratio profile R(A 3 ,z) and its 
uncertainty 6 R(A 3 ,z). The estimate for R(A 3 ,z) is given by 

-e(A ) e(X 3 ) 

R(A 3 ,z) = 1 + [R(Xi,z) - l] v F(X 3 ,Ax,z)A 1 X 3 > (x 

where f is a best-estimate for the ratio of particulate back- 

scattering at wavelength A q to that at A,, and e(A) is the 

_ 8,9 

wavelength exponent for gaseous elastic backscattering. 

[ 6 R(A 3 ,z) is defined in Section III-A. ] 

Step _3: Use the estimated short-wavelength scattering ratio profile 

RQ^z) and the measured short-wavelength signal profile 
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A 1 = 1.06 fj.m \ 2 ~ 0-53 firri \ 3 = 0.35 /im 



FIGURE 1 MULTIWAVELENGTH ANALYSIS PROCEDURE TO RETRIEVE PROFILES 
OF SCATTERING RATIO AND MOLECULAR DENSITY 

Numbers indicate steps described in the text. 


S(A 3 ,z) to obtain the lidar-derived density profile and its 
uncertainty. The relative density profile is 

2 

(z - z L ) S(A 3 , z ) 

D(z) 5 _ , 

Q ( a 3» 2 » 2 l) r C a 3» z ) (2) 

which follows from Eq. (2) of Part Two of this report, where 
the symbols are defined. If absolute density p is known at 
any height z, the absolute density profile can be obtained as 

D(z) 

pL^ 2 ) = p( 2 ) . 

D(z) (3) 
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In this paper we usually choose z to be an assumed isopycnic 
level, where the density can be estimated to high accuracy ji 
priori , because of the very small natural variability there* 
The uncertainty in Eqs. (2) and (3), including the uncertainty 
in assuming an Isopycnic level, is discussed in Section III-A. 

Step ki Substitute the lidar-derived density p^(z) for the conven- 
tional density p c (z) at all heights z where 6p^(z) < <$p c (z). 
Retain the conventional values elsewhere. Call the resulting 
profile the "composite" density profile. 

Step _5: Repeat Step 1 using the composite density profile and Its 

uncertainty in place of the conventional density profile and 
its uncertainty. Also, if other lidar wavelengths are avail- 
able (e.g., 0.53 vim in Figure 1), use the composite density 
profile with the signal profile at each other wavelength to 
derive the corresponding scattering ratio profiles. 

o 

The key to this rather simple procedure is that the ratio of parti- 
cle to gas backscattering depends very strongly on wavelength. For 
example, the approximate wavelength dependence of gas backscattering 

yields 

B g (0 .355 ym) /B g ( 1 .064 pm) * 81, (4) 

whereas a typical wavelength dependence of particle backscattering, 
with typically* * » 


0 < ct < 2 , 


yields 


Bp(0.355 p.m)/Bp( 1 .064 pm) a 1 to 9. 


(5) 


o 
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This combination yields backscatter mixing ratios B p /B g that are much 
larger at the longest wavelength than at the shortest wavelength. For 
example, combining Eqs. (4) and (5) yields 

B p (0.355 pm) B p (1.064 ym) 

» (0.012 to 0.11) x . 

Bg(0.355 ym) Bg(1.064 ym) (6) 

Figures 2 through 4 of Part Two of this report show that, below 40 km, 
Eq. (6) applies for all aerosol and cloud components of each model atmo- 
sphere. This is significant because the aerosol backscatter mixing 
ratios shown in that figure were derived from Mie-scattering calcula- 
tions for the size distributions and refractive indices appropriate to 
each model component (mesospheric, stratospheric, upper tropospheric, 
Saharan, and marine) rather than for an assumed wavelength dependence, 
such as the used above for illustration. (Unfortunately, we know of 
no measurements of natural aerosol and cloud particulate backscatter 
spanning the wavelength range 0.35 to 1.06 ym.) 

This typical wavelength dependence of backscatter mixing ratio 
means that an aerosol or (subvisible) cloud layer that yields a 0.355-ym 
scattering ratio of 1.01 (hence perturbing the lidar signal profile by 1 
percent) typically yields a 1.064-ym scattering ratio of 1.09 to 1.81; 
hence, this layer is easily detectable in the 1.064-ym profile, even if 
errors of 2 to 5 percent are present in the conventional density profile 
used in Step 1. Detection of such a particulate layer alerts the user 
to the possible contamination of the lidar-derived density profile and, 
through Step 2, allows for an approximate, yet useful, correction for 
these effects. Section III quantifies these arguments and shows how the 
indication of particulate contamination is automatically included in the 
error bar <5p L . 

Some possible refinements in the multiwavelength solution scheme of 
Figure 1 are immediately apparent. For example, one could iterate, fol- 
lowing Step 5 with several repeats of Steps 2-5, or one could include 
information on R(A 2 - 0«53 ym) in Step 2 when estimating R(X^) from 
R(Xj). However, the stability of such an iterative scheme is not 
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guaranteed, and the larger relative errors in B p (A 2 ) /B g (X 2 ) [see Part 
Two of this report] may actually degrade the accuracy of estimated 
R(A^), unless care is taken to exclude results from X 2 that do not 
satisfy some error criterion. For these reasons, and also because our 
simulations indicate that the basic scheme in Figure 1 can yield very 
useful results in many situations, we have limited the analyses in this 
paper to the basic scheme. Refinements of the basic scheme are probably 
best left until new insights have been developed by applying the basic 
scheme to actual measurements. 


B. Pressure and Temperature Profiles 

To derive pressure and temperature profiles from the density pro- 
file obtained above, we combine the hydrostatic relation 

dP(z) = -g(z)p(z)dz (7) 


with the ideal gas law 

P(z) 

T(z) - T' v (z) = . 

Cgp(z) (8) 


In these equations P is pressure, g is the acceleration of gravity, 
is the gas constant for dry air, T is the absolute air temperature, and 

1 O 

T' v is the adjusted virtual temperature. Differences between T' v and 
the actual temperature T are caused by the presence of water vapor; how- 
ever, at heights above ~5 km (to which we will find our density mea- 
surement technique is limited in practice), typical relative humidities 

1 3 

and temperatures yield 

It - T' v | < 1 K, (9) 


125 


I 


and so we henceforth approximate the two as being equal. Since Eq. (7) 
is in practice applied to the finite interval 

Azi,i±i s Z1 - z i±1 (10) 


where i identifies a lidar range bin, and since p(z) is increasing or 
decreasing exponentially in this interval, we explicitly integrate Eq. 
(7) over the finite interval by assuming a constant scale height 

in this interval. This yields 


APi,i±l 



g(z)|o(z i )e 


-(z-z i )/H 1>i±1 


dz 


(ID 


= g( ZjL ,z i±1 ) H i>i±1 [p(z i± i) - p(zi)] , 

where is the scale height, obtained from 

H i,i±l = < z i±l - z i )/ln[p(z i )/p(z i±1 )] . 

Note that Eq. (11) is equivalent to 

APi,i±l = -g(zi,z i± i)p(zi,z i± i) Azi,iii , 

provided we define 

p(z i± i) - p(zi) 

p(Z£,Z-^j-2) =: “ ' . 

ln[p(zi)/p(zi ± i)] 

Use of Eq. (15) in place of the linear average 

p(zi,zi ± i) = 0.5 [p(zi) + p(z i±1 )] 


( 12 ) 


(13) 


(14) 


(15) 


(16) 
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improves accuracy significantly when Az^ is an appreciable fraction 
of the scale height The choice of a linear or logarithmic 

definition for the interval-average gravitational acceleration g is not 
critical, because g(z) is a very slowly varying function of z. 

If we now define a reference height 2 at which the pressure or tem- 
perature can be guessed with an accuracy 6P or 6 t, we have 

P(z x ) = F + APi , (17) 


where 


A?l = “g(zi,2)|5(zi,Z)AZi 


(18) 


As I = zj - z 


(19) 


and I, 1+1 are the Indices of the lidar data heights surrounding 2. The 
complete pressure profile can then be generated using 


i 

p(z ± ) = P( 2l ) + £ APj,j-i 

j-i+i 


( 20 ) 


for i > I and 


P(zf) = P(z X ) + 2 APj, j+1 

j=I-l 


( 21 ) 


for i < I. Substituting Eqs. (14), (15), and (17) through (19) into Eq. 
(20) yields, for the complete pressure profile, 


P(zi) = P + Kp(2,z X ) + XJ k j»j~l p( z j» z j-l) 

j«I+l 


( 22 ) 
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with an analogous expression for i < I. Here, we have defined 


E = -g(zi,z)A 2 I » 


(23) 


k j,j±l - s(zj, z j±l)Azj, j±l 


(24) 


The temperature profile is obtained by substituting Eq. (22) into 
Eq. (7) to yield 

P E p(z,zj) * P^ z j ,z j~l^ 

T(zi) = + + 2 • 

c g p( z i) C g p( Zi ) j=I+1 C g p( Zi ) (25) 


If a reference temperature T, rather than a pressure P, is specified, it 
is convenient to rewrite Eq. (25) as 

p E p(2, Z j) * p( Z j, Z j-i) 

T( Zi ) = T + + X . 

P (z i> C g p( Zi ) j„ I+1 C g p( 2l ) (26) 


where p is obtained by exponential interpolation between pi and px+i« 
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Ill ERROR ANALYSIS 


A. Gas Density and Particle Backscatter Profiles 


Applying standard error-propagation techniques^ to Eq. (1) yields 


'6r 3 N 


r 3 


'r 3 - 1 N 

r 3 


<5Ri 


'5^ 


3,1 


R 1 " 1 


^3, 1 


(27) 


where we have introduced the shorthand notations 


Rn = R(A n ,z) , 


(28) 


^m,n = ^(^m*^n> z ) • 


(29) 


Similarly, one obtains from Eqs. (2) and (3) 

2 2 2 


PL 


PL 


¥ 


+ 


S(D/D) 


D/D 


(30) 


where 


2 2 2 2 2 2 
'<5(D/f>)\ / 6 s 3\ /^® 3 \ /5Q 2 (X 3 ,z,z)\ /6R 3 \ /<$R 3 ) 


D/D 


S 3 


s 3 


Q (A 3 , z , z) 


r 3 


R 3 / (31) 


with the shorthand notations 


P " p(z> , 


(32) 
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S m S(A. m> z) , Sjq S(A m> z) 


(33) 




(34) 


Also, we have made use of the fact that 


2 2 2 
Q (A,z,z l )/Q (A,z,z l ) = Q (A,z,z) = exp 


-2 J E(X ,z' )dz / 

z 


(35) 


[See Eq. (8) of Part Two and discussion.] 

Equation (31) neglects a negative term proportional to the covari- 
ance of R 3 and R^. This covariance is nonzero because R 3 and R 3 are 
positively correlated through their dependence on and also because 

Rl and R^ are positively correlated through their dependence on 
P^in^i), p c (z*), and S(X,z*) [cf. Eq. (4) of Part Two and Eq. (27) of 
this part.] However, an algebraic solution for the covariance term is 
extremely complicated, and its magnitude turns out to be small for the 
practical cases we will simulate. Because the neglected covariance term 
in Eq. (31) is negative, Eqs. (30) and (31) give a slight overestimate 
for and are thus conservative. The simulations in Section IV 

shed light on the magnitude of the neglected covariance for practical 
cases. 


B. Temperature Profiles 

In assessing uncertainties in the lidar-derived temperature pro- 
file, we make use of the fact that density p occurs only as a ratio in 
Eq. (26) and in all terms but the first in Eq. (25). As a result, 
lidar-derived T is independent of any constant multiplier in the lidar- 
derived density profile — particularly of the lidar-derived factor D and 
the assumed isopycnic density (6 [see Eqs. (2) and (3)]. [This is true 
at all heights if the parameter specified at the reference height z is 
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temperature; if instead the reference parameter is pressure, T depends 
on p/D near and above z because of the first term in Eq. (25).] 


Substituting Eqs. (3) and (2) into Eq. (26) and approximating den- 
sity averages by the linear form in Eq. (16) yields 

2 

/D^ 



1 1 

1 

2 C g | 

/ 

Di 


K(D' + D x ) 


i-1 


+ E 

j=i+l 


k j> j-l 


(D j-l + Dj) + kiji-i T>±-l 


+ k 


i,i-l 
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When pressure P is specified at the reference height, we still can 
derive Eq. (36) from Eq. (25), except that the first term is replaced by 

P t) 

C g Di p (38) 

The primed relative densities Dj defined by Eq. (37) are identical to 

the unprimed quantities defined by Eq. (2), except that all two-way 
2 

transmissions Q are referred to the height z^ where temperature is 
being computed. This definition removes the effect of errors in 
transmission over all heights except those between z^ and zj. (Recall 
that all Zj are between and 2.) This removal occurs because 
transmission errors at all other heights affect D^ and Dj equally, and 
thus cancel in the density ratios of Eqs. (25) and (26). Hence, 
although D^ and Dj are correlated in this respect, D^ and Dj are not. 
Writing Eq. (36) in terms of primed relative densities therefore 
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simplifies the error analysis by replacing a correlated set of variables 
by a set that has this correlation removed. 

14 

Applying error propagation analysis to Eqs. (36) and (38) yields 
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If the parameter at the reference height z is temperature, then 
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if instead pressure is specified, then 
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IV SIMULATION PROCEDURE 


To test the algebraic error expressions described above, and also 
to develop insights and data reduction software, we extended the simula- 
tion procedure described in Part Two. The extended procedure follows 
the analysis flow shown in Figure 1. In Steps 1 and 5, which are con- 
ventional single-wavelength retrievals, random errors are introduced 
into the simulated measurements and retrievals as described in Part 
Two. (See Figure 1 of Part Two.) In Step 2, R(A^,z) is computed from 
Eq. (1). Currently we are using 

-1 

B p ( x 3 ,z) / A 3 \ 

^ ( A 3 , Aj ,z) = — I j , 

B p (>i,z) \A]/ (44) 

independent of height, in this step. The "random errors" in this step 
are provided by the differences between Eq. (44) and the wavelength 
dependence of the input aerosol and cloud models at each height. 

[Recall that multiwavelength backscattering for each aerosol submodel is 
calculated directly from size distribution and refractive index, rather 
than from an assumption, such as Eq. (44); the resulting height and 
latitude dependence of ¥ can be seen in Figures 2 through 4 of Part 
Two. ] 

In Step 3 [Eqs. (2) and (3)], random errors in signal and transmis- 
sion are introduced, just as in the single-wavelength simulations 
described in Part Two. In addition, a random error is added to the 

a 

model value for p(z) to simulate an incorrect guess for density at the 
assumed isopycnic level. 

No random errors are added in Step 4. 
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In Step 5, care is taken to retain the same errors in the conven- 
tional density profile as were used in Step 1, as well as the same sig- 
nal, transmission, and normalization errors for A^. 



V SIMULATION INPUTS 


The simulations shown in this paper use the same atmospheric 
models, lidar parameters, background lighting, and error sources for 
signal, conventional molecular density, transmission, and R m ^ n as used 
in Part Two, In addition, we have used 

<5 y 

= 100 % 

4/ (45) 

and 

( S P\ 

I 1 = 2%, z = 8 km 

V P / (46) 

as measures of the lcr uncertainty in the assumption of Eq . (44) and in 
estimating an ,f isopycnic M density, respectively. Equation (45) was 
obtained by inspecting the range of ¥ values, and their rms deviation 
from Eq. (44), for the aerosol and cloud models described in Appendix A 
of Part Two. Equation (46) is a very rough estimate obtained by 
inspecting a collection of model molecular density profiles, plus a 
small set of balloon-measured temperature, pressure, and density pro- 
files. We recommend that Eq. (46) be reevaluated using a large set of 
radiosonde data. However, for the present we note that retrieved 
scattering ratio profiles and temperature profiles are independent of p, 
so Eq. (46) has no effect on these products or their uncertainties. 

[Note the discussion just before Eq. (36) for the exception to this case 
for temperature profiles.] However, retrieved particulate backscattering 
profiles, density profiles, and pressure profiles do depend on p, and 
hence their derived uncertainties are sensitive to Eq. (46). 
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VI SIMULATED PERFORMANCE 


A, Tropical , Nonvolcanic , Cloud-Free , Saharan ; 

Nighttime ; Ax ** 200 km 

Figure 5(a) of Part Two^ has shown the profile of 1.06-ym back- 
scatter mixing ratio (Bp/Bg = R - 1) obtained from a simulated nighttime 
retrieval that combined the 1.06-ym signal profile with a conventional 
density profile- That retrieval is Step 1 of the multiwavelength 
analysis procedure (see Figure 1 of this paper). Step 2 uses the 1.06- 
|jm scattering ratio profile (and its uncertainty) with Eqs. (1), (27), 
(44), and (45) to estimate a 0.355-ym profile of scattering ratio (or 
backscatter mixing ratio) and its uncertainty. The result is shown in 
Figure 2(a). The estimated profile (dots) is directly proportional to 
the 1.06-yra profile in Figure 5(a) of Part Two, because of the assump- 
tion of a height-independent 'P in Eq. (44). Also, the error bar on each 
dot in Figure 2(a) extends at least from zero to twice the dot, because 
of the assumed ±100 percent error in H'. [See Eqs. (45) and (27).] 

Comparison of the estimate (dots) with the model (solid line) in 
Figure 2(a) shows that for the aerosol model in this scenario, Eq. (44) 
leads to a systematic underestimation of R(0.355 ym) in the stratosphere 
and a systematic overestimation in the lower troposphere. This suggests 
that errors could be reduced by using a height-dependent 4* that assumes 
a steeper wavelength dependence for particulate backscattering in the 
stratosphere, and a less steep one in the troposphere. However, con- 
firming this would require careful inspection of large sets of size 
distribution and composition data, plus associated Mie-scattering calcu- 
lations, or, better yet, analysis of a representative set of mul- 
tiwavelength backscatter measurements. This is beyond the scope of the 
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FIGURE 2 LOW-LATITUDE NIGHTTIME SIMULATION RESULTS, 0.355 /im 

(a) Backscatter mixing ratio profile inferred from 1. 064-Aim profile in Figure 5(a) 
of Part Two; (b) Signal and background profiles; (c) Density profile retrieved from 
(a) and (b), expressed as ratio to model; (d) Temperature profile derived from (c), 
compared to model. 
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current study; moreover, the simple forms in Eqs. (44) and (45) provide 
quite useful results and also demonstrate the analysis method. 

Step 3 (see Figure 1) combines the estimated backscatter mixing 
ratio profile in Figure 2(a) with the measured signal profile [shown for 
a single shot in Figure 2(b)] to obtain the lidar-derived density pro- 
file given by Eqs. (2) and (3). The simulated result is shown in Figure 
2(c), expressed as a ratio to the exact (model) density profile. Notice 
that the lidar-derived profile is biased to the right. This results 
from a 2 percent error in the value guessed for p at the assumed iso- 
pycnic level (8 km), and a partially offsetting error in signal S at the 
same height. [See Eqs. (2), (3), and (30) through (35).] The error bars 
in Figure 2(c) contain only the uncertainty in the relative profile D(z) 
as given by Eq. (31). Both the error bars and the scatter among the 
dots show that the lidar-derived relative density profile for this simu- 
lation has a relative lo- error of ~1 percent or less for heights between 
8 and 20 km. This is considerably better than the 2 percent assumed for 
the conventional density profile and, moreover, it comes in a height 
region where reduction of density errors is most strongly needed to 
improve particulate retrievals (see below). 

Figure 2(d) shows the temperature profile derived from the relative 
density profile in Figure 2(c) using Eq. (26) and an estimated reference 
temperature T at 42 km. In the case shown, T was in error by +8 K, and 
also D(z) was in error by ~+2 percent and ~-3 percent at 37 and 33 km, 
respectively (because of lidar signal measurement errors). However, 
errors in retrieved temperature rapidly become independent of these 
high-altitude errors, and the errors decrease to between ±1 and ±2 K 
between 20 and 8 km. [This can be most easily seen in the expanded 
scale on the right of Figure 2(d).] These small errors, combined with 
lidar 's excellent vertical resolution (0.5 km in this height region and 
scenario), give an excellent retrieval of the tropopause structure. 
Existing spaceborne nadir-viewing sensors are unable to perform such 
tropopause retrievals. 
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Effects of the tropical nonvolcanic stratospheric aerosol layer can 
barely be seen in Figure 2(c), especially between 22 and 27 km, where 
the underestimate in R(0.355 ym) is greatest [see Figure 2(a)], How- 
ever, these effects are comparable in magnitude (~1 percent) to the 
signal-measurement errors, and so are hard to distinguish from the sig- 
nal error. Effects of the strong Sahara and marine aerosol layer are 
very obvious below 7 km, leading to extremely large errors in the 
lidar-derived density and temperature profiles. However, the error 
bars, calculated from Eqs. (30) and (39), include the effects of these 
aerosol layers (detected in the Step 1 analysis at 1.06 ym) and immedi- 
ately warn of their presence. 

Figure 3 shows results of aerosol retrievals at 1.064 and 0.532 ym 
using the lidar-derived density profile [Figure 2(c)] above 7 km and the 
conventional density profile below. These results can be compared with 
Figure 5 of Part Two, which was obtained from the same signal mea- 
surements but used the conventional density data at all heights. (Simu- 
lated conventional density rms errors were 2 percent below 30 km and 3 
percent above.) This comparison shows that the 0.53 ym retrieval that 
uses the lidar density profile [Figure 3(b)] is considerably better in 
the upper troposphere and lower stratosphere (~8 to 20 km), where den- 
sity measurement errors have their largest effect (because of the small 
backscatter mixing ratios). Despite these small backscatter mixing 
ratios, this is an important region for aerosol measurements, because 
absolute backscatter and extinction here are comparable to those in the 
stratospheric peak. Hence, this region contributes significantly to the 
total stratospheric/upper tropospheric optical depth, and erroneous mea- 
surements here can significantly degrade the accuracy of derived optical 
thickness and column backscatter values. Studies of stratospheric- 
tropospheric exchange also require accurate measurements in this region. 
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FIGURE 3 LOW-LATITUDE NIGHTTIME SIMULATION RESULTS, 1.064 AND 0.532 Mm, 

USING LIDAR DENSITY DATA ABOVE 7 km, CONVENTIONAL DENSITY BELOW 

(a,b) Backscatter mixing ratio profiles; (c) Relative uncertainty in particulate 
backscattering broken down by source. 


141 





B 


Midlatitude , Volcanic , Cloudy , Marine ; Nighttime ; ^x * 200 km 


Figure 4 of this paper shows 0.355-ym results for the same scenario 
as Figure 7 of Part Two. These results extend only to 44 km because 
signal-measurement errors exceed 3 percent above that height. In this 
case the lidar density profile [Figure 4(c)] was again normalized to an 
assumed isopycnic layer at 8 km, with an offset error of +2 percent. 
However, in this case the offset increases with increasing height 
because of transmission errors that accumulate through the considerable 
optical thickness of the upper two cloud layers and the moderate vol- 
canic stratospheric aerosol layer. This leads to a bias error of about 
5 percent above 20 km. Because this error is larger than that expected 
for the conventional density profile (assumed 3 percent for this case, 
to simulate use of a local seasonal model density profile), a more accu- 
rate absolute lidar density profile could have been obtained by normal- 
izing to the conventional density somewhere between 20 and 30 km (where 
lidar signal errors are less than 2 percent). However, we have retained 
the normalization at 8 km to demonstrate that this type of bias or 
offset error does not significantly affect temperature profile 
retrievals above the layers in which the transmission errors occur. 

[Note that the density error bars in Figure 4(c) do not include the den- 
sity offset error, and thus are not expected to span the difference 
between the dots and the vertical line. They do describe well the rms 
difference between each dot and the normalization dot at 8 km, as 
intended. ] 

Figure 4(d) shows the temperature profile obtained by integrating 
the lidar density profile downward from 44 km, where the reference tem- 
perature was guessed with an error of +8 K. As in Figure 2(d), after 
integrating downward by _8 km to become independent of the reference 
temperature guess and the larger signal-induced density errors, the tem- 
perature profile is retrieved to an accuracy of 1 to 3 K (between 32 and 
15 km), in spite of the ~5 percent density offset error above 20 km. 

The reason for this was explained between Eqs. (38) and (39). Physi- 
cally, the reason is that pressures derived from the lidar density 
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FIGURE 4 MID-LATITUDE NIGHTTIME SIMULATION RESULTS, 0.355 jum 

(a) Backscatter mixing ratio profile inferred from 1.064 jLtm profile in Figure 7(a) 
of Part Two; (b) Signal and background profiles; (c) Density profile inferred from 
(a) and (b), expressed as ratio to model; (d) Temperature profile derived from (c), 
compared to model. 
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profile are about 5 percent too large between 32 and 20 km; then, in 
using the gas law [Eq. (8)] to obtain temperature, these pressures are 
divided by a density that is also too large by about 5 percent. 

Between 21 and 18 km, the moderate volcanic stratospheric aerosol 
introduces a local error of _+3 percent in density and -3 K in tempera- 
ture. (The maximum relative temperature error, -3/215 = -1.5%, is less 
than the maximum relative density error because of a partial compensa- 
tion in pressure error caused by aerosols just above the peak.) Below 
the stratospheric aerosol peak, positive temperature errors occur as the 
transmission-induced density bias errors decrease and there is a 
mismatch of pressure and density errors. Below 12 km, the subvisible 
and visible cloud layers cause very large temperature and density 
errors . 

For this scenario, use of the lidar density profile above 12 km 
produced only minor improvement in backscatter mixing ratio retrievals 
at 0.53 and 1.06 pm (not shown). This is a combined result of the rela- 
tively large lidar density errors and the relative insensitivity of the 
conventional retrievals [Figure 7 of Part Two] to density errors, 
because of the relatively large moderate volcanic backscatter mixing 
ratios. This is an example of the general rule that, when particulate 
layers significantly perturb the lidar density profile, high-accuracy 
density values are not required to obtain accurate particle retrievals. 
On the other hand, when backscatter mixing ratios are small (thus 
requiring highly accurate relative density profiles for particle 
retrievals) the lidar density profile has very small particle-induced 
errors. 

C. High-Latitude , Nonvolcanic , Cloud-Free , Marine , Nighttime , 

/\x » 2000 km 

Figure 5 shows 0.355 ym results for the same scenario as for Figure 
9 of Part Two. In this case, results extend up to 56 km, because the 
longer integration time (280 s) reduces signal measurement errors to 
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FIGURE 5 HIGH-LATITUDE NIGHTTIME SIMULATION RESULTS, 0.355 Atm 

(a) Backscatter mixing ratio profile inferred from 1.064-jum profile in Figure 9(a) 
of Part Two; (b) Signal and background profiles; (c) Density profile inferred from 
(a) and (b), expressed as ratio to model; (d) Temperature profile derived from (c), 
compared to model. 
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about 3 percent at that level. Again, the lidar density profile was 
normalized at 8 km with a bias error of +2 percent; however, most of the 
profile has little offset because of a fortuitously large positive sig- 
nal error at that height. (As shown above, results for temperature and 
aerosol retrievals would have been very similar without this fortuitous 
signal error, because of the insensitivity of these results to bias 
errors.) The relative density profile has rms errors of 0.5 to 2 per- 
cent between 4 and 56 km. These errors are caused primarily by signal 
errors, although the effect of the nonvolcanic stratospheric aerosol 
peak can barely be seen at 15 km. 

The lidar-derived temperature profile [Figure 5(d)] was obtained by 
integrating down from 56 km (i.e., near the stratopause) , where the 
reference temperature again had an error of +8 K. As before, tempera- 
ture errors had decreased to „3 K at a height (47 km) that was 9 km 
below the reference height. Between 4 and 45 km, rms temperature errors 
are 1.2 to 2.5 K, with the exception of a -3.8 K error at 15 km caused 
by the nonvolcanic stratospheric aerosol layer. 

Figure 6 shows 1.06- and 0.53-ym aerosol retrievals using the lidar 
density profile above 5 km and the conventional profile below. Com- 
parison with Figure 9 of Part Two shows considerable improvement in the 
ability of the 0.53-ym retrieval to replicate the vertical structure of 
the stratospheric aerosol. The 1.06-ym retrieval also improves signifi- 
cantly. 
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FIGURE 6 HIGH-LATITUDE NIGHTTIME SIMULATION RESULTS, 1.064 AND 0.532 fim, 

USING LIDAR DENSITY DATA ABOVE 5 km, CONVENTIONAL DENSITY BELOW 

(a,b) Backscatter mixing ratio profiles; (c,d) Relative uncertainty in particulate 
backscattering broken down by source. 
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VII SUMMARY AND CONCLUSIONS 


We have described an analysis technique that combines lidar signal 
profiles at long and short wavelengths (in this case, 1.064 and 0.355 
pm) to yield profiles of density and temperature, and aerosol and cloud 
backseat ter, each with an associated uncertainty profile. The technique 
begins with conventional analysis of the long-wavelength profile, which 
leads to an estimate of particle backscatter at the short wavelength. 

This estimate is combined with the short-wavelength signal profile to 
yield a relative density profile and an uncertainty profile that 
includes effects of aerosol contamination. 

This relative density profile is converted to an absolute density 
profile by normalization at an assumed isopycnic level or a level where 
conventional density data are available to good accuracy. The lidar- 
derived density profile is used to improve aerosol and cloud retrievals 
at 0.53 and 1.06 pm, and is also integrated vertically (from a reference 
height where temperature or pressure are guessed) to yield pressure and 
temperature profiles. The two keys to the success of this analysis pro- 
cedure are: 

• Returns at 0.355-pm are much more sensitive to gas backscatter- 
ing than are those at 1.06 pm. (Hence, on a relative basis, 

0.355 pm returns are much less sensitive to aerosol contamina- 
tion, whereas 1.06 pm returns readily indicate regions where 
particulate contamination might be significant.) 

• Aerosol retrievals and the lidar-inferred temperature profile 
are relatively insensitive to bias errors in the density pro- 
file, so that errors in absolute density normalization have lit- 
tle effect on these data products. Also, this reduces the 
effect of some transmission errors. 

Error analysis equations were developed to quantify the above argu- 
ments and to provide error bars on all retrieved quantities. Also, to 
check the algebraic error expressions, the numerical simulation code was 
extended to inject appropriate errors into the multiwavelength retrieval 
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procedure. In addition to the four error sources considered for 
single-wavelength, conventional retrievals (i.e., errors in signal, 
transmission, lidar calibration, and conventional density profiles), we 
also included errors in absolute density normalization, in the reference 
temperature or pressure, and in estimating short-wave length particle 
backscatter from a long-wavelength measurement. 

Simulations were run for the low, middle, and high-latitude model 
atmospheres. In general, cirrus clouds (subvisible as well as visible) 
and lower tropospheric aerosols introduced larger errors in density and 
temperature than those in conventional data; daytime signal measurement 
errors were similarly large in the stratosphere and above. Hence, use- 
ful lidar density and temperature measurements were restricted to night- 
time conditions in the stratosphere and (cloud-free) upper troposphere. 
(Signal errors also prevented useful measurements in the mesosphere if 
horizontal resolution was 2000 km or less.) 

Simulations were run with vertical resolutions of 0.5 to 1.0 km in 
the upper troposphere and lower stratosphere, and 2.0 km in the upper 
stratosphere. Horizontal resolution of 200 km provided useful signal 
accuracies (0.5 to 2.0 percent) below 40 km; however, this had to be 
increased to 2000 km to extend useful signals to 55 km. Retrieved rela- 
tive density profiles had rms errors of 0.5 to 2 percent in the upper 
troposphere and stratosphere. The limiting error at most heights was 
signal measurement error. However, nonvolcanic stratospheric aerosols 
introduced errors of 1 to 2 percent at the mixing ratio peak (typically 
~27 km at low latitudes to ~16 km at high latitudes). For a given 
amount of nonvolcanic stratospheric particle backscatter, the density 
error tends to be larger at low latitudes, because the low-latitude peak 
occurs at greater heights, where the smaller absolute density yields 
larger relative errors. A midlatitude, moderate volcanic stratospheric 
aerosol model introduced relative density errors of about 3 percent at 
the peak (-18 km) • 
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Use of the lidar-derived relative density profiles in place of con- 
ventional profiles significantly improved the accuracy of aerosol 
retrievals in the nonvolcanic stratosphere and upper troposphere — 
especially at 0.53 pm, where conventional density errors gave poor 
retrievals at these heights. The impact of lidar-derived density errors 
on volcanic stratospheric aerosol retrievals was minimal, both because 
of the insensitivity of such retrievals to density errors and because of 
the poorer accuracy of the volcanic lidar density profile. 

For cloud-free, nonvolcanic conditions, lidar-derived temperature 
profiles had ms errors of 1.2 to 2.5 K in a layer bounded on the bottom 
by strongly scattering tropospheric aerosol or cloud layers and on the 
top by a height -.8 km below the reference height where temperature or 
pressure was guessed. For the cloud-free model atmospheres used here, 
this lower bound was ~5 km; the upper bound was ~32 km for a horizontal 
resolution of 200 km and ~47 km for a horizontal resolution of 2000 km 
(i.e., the reference heights were at _40 and ~55 km, respectively). The 
midlatitude, moderate volcanic stratospheric aerosol introduced tempera- 
ture errors of ~3 K between 17 and 20 km. These simulations assumed an 
8 K rms error in the temperature guessed at the reference height (40 to 
55 km) . 

The major advantage of these lidar-derived temperature profiles 
over those obtained by passive nadir-viewing spaceborne sensors is their 
fine vertical resolution (—0.5 km in the upper troposphere and tropo- 
pause region, increasing to Jl km in the upper stratosphere). This 
resolution would permit us to define the tropopause and temperature-wave 
structures to a degree of detail never before possible from space. 

In general, the simulations validated the algebraic expressions 
used to put error bars on lidar-derived density and temperature pro- 
files. This is important because the errors reported above apply to 
specific atmospheric situations and will, for example, increase consid- 
erably whenever strong particulate (cloud or aerosol) layers are encoun- 
tered. In this respect the outstanding feature of the multiwavelength 
analysis technique is that it simultaneously retrieves density and 
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particle profiles, and includes particle contamination effects in the 
density and temperature error bars. Thus, when significant cloud or 
aerosol contamination occurs, the user is immediately aware of it. 
Overcoming these strong particle contamination effects requires more 
sophisticated techniques, such as the high spectral resolution method of 
resolving the narrow particle-backscatter line from the Doppler- 
broadened Rayliegh line,^*^ or differential-absorption measurements of 
a well-mixed gas, such as molecular oxygen. 

1 8 

A previous contract report also describes the two-wavelength 
elastic backscatter technique for particle/gas separation as a means of 
measuring gas density profiles. Although that report is based on the 
same ideas as used here, it differs in the following ways: 

• The error-propagation equation (p. 101, Ref. 18) has a fundamen- 
tally different form from our equations and does not include any 
calibration error terms [the SR^n* ^S*, 5D*, SR, SS, 5p of our 
Eqs. (7), (9), (11), (12), Part Two, and (30) and (31), Part 
Three]. The Ref. 18 equation evidently results from assuming a 
perfectly calibrated lidar, which does not need to use a ’’clean" 
or "isopycnic" atmospheric layer to relate detector output to 
backscatter coefficient or molecular density. Our experience 
with actual lidar backscatter measurements has been that some 
sort of atmospheric calibration is almost always required. This 
would be especially so for a spaceborne system, where a cali- 
brated target would not be available to measure any system 
degradation that occurred in orbit. Whereas the error- 
propagation equation of Ref. 18 shows many of the relevant error 
processes, we do not feel that it is directly applicable to 
practical lidar measurements. 

• The only signal-measurement error considered is the JW pulse 
counting error. Other possible signal-measurement errors (e.g., 
current saturation, nonlinearity, pulse-to-current splicing), 
which led to the C 6S/ S ) m -^ T1 term in this study, are neglected. 

• No simulations with actual lidar data-processing algorithms and 
random number generators that simulate each error source are run 
to check the algebraic error-propagation equation. 

• The derived expected molecular density error between 5 and 30 km 
(p. 109 of Ref. 18) varies from 4 to 14 percent, with a peak at 
~18 km, provided one assumes an uncertainty in the wavelength- 
dependence ratio ^ [= B (0.355 ym)/B (1.064 ym)] that we con- 
sider reasonable [i.e., 50 to 100 percent — see Eq. (45) and Sec- 
tion III-B of Part One]. These expected errors considerably 
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exceed those (0,5 to 2 percent) obtained in this study for non- 
volcanic model atmospheres, even though this study assumed = 
100 percent. 

The explanation for the final discrepancy is the different aerosol 
models used in the two studies. Reference 18 used the LOWTRAN 3B aero- 
sol model. ^ In turn, LOWTRAN 3B uses the Elterman 1968^ vertical pro- 
file of midvisible particle extinction in the upper troposphere and 
stratosphere. As shown by Shettle and Fenn, the Elterman 1968 extinc- 
tion profile slightly exceeds the moderate volcanic extinction profiles 
of Shettle and Fenn. In addition, Ref. 18 assumes a backscatter phase 
function of 0.04 sr ^ for all wavelengths. As shown by Figure B-l of 
Part Two, this value may be appropriate for some tropospheric aerosols, 
but it is a factor of two or more too large for background stratospheric 


aerosols. The 

result is that 

the Ref. 18 

stratospheric 

particle 

back 1 

scatter values 

significantly 

exceed our moderate volcanic values. 

and 

they exceed our 

nonvolcanic values by factors between 6 

and 9 at 


wavelength 0.353 pm. The model values of 

B p (X)/B g (x) at 

the statos- 

pheric peak may 

be summarized 

as : 



Wavelength 

(ym) 

High- 

Latitude 

Nonvolcanic 

Mid- 

Low- Latitude- 

Latitude Moderate 

Nonvolcanic Volcanic 

Reference 

18 

0.355 

0.02 

0.03 

0.09 

0.17 


0.532 

0.05 

0.1 

0.25 

0.45 


0.694 

0.1 

0.25 

0.7 

0.85 


1.064 

0.3 

1.0 

2.0 

2.4 



The middle three columns are taken from Figures 2 through 4 of Part Two. 
The final column is from pp. 97 and 99 of Ref. 18, with the 0.532 and 
0.694 ym wavelength values logarithmically interpolated between the 
0.355 and 1.064 pm values. 
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This comparison raises the questions of 

(1) How well are the various models supported by measurements? 

(2) How often do "nonvolcanic" and "moderate volcanic" conditions 
occur? 

Our models are based on dustsonde measurements of particle number, cou- 
pled with model size distribution and refractive index ranges that are 
also based on measurements. Our models have been validated by comparis- 
ons to ruby lidar backscatter measurements (X = 0.694 ym) , but when the 

models were developed no 1.064- or 0.353-ym backscatter measurements 

22 

were available for validation. Since then, Shibata et al. have pub- 
lished nonvolcanic profiles of B /B at 1.06 ym. Although 0.353-ym mea- 

P 8 

surements were not made, 0.532-ym measurements were made, and these are 
useful in estimating upper bounds. Ruby and Nd-YAG^ measurements 

of peak stratospheric B (x)/B (A) may be summarized as: 

r o 


Wavelength 

(ym) 

High- 

Latitude 

Nonvolcanic 

Low- 

Latitude 

Nonvolcanic 

Mid- 

La titude- 
Moderate 
Volcanic 

0.532 


0.08-0. I 22 


0.694 

0.1 23 

0.2 24 

0.7 25 ’ 26 

1.064 


0.25-0.4 22 



These measured values agree quite well with the corresponding model 
values summarized above, but they are markedly less than the Ref. 18 
model values. 

Regarding the question of frequency of occurrence of "nonvolcanic" 

and "moderate volcanic" conditions, the time series of ruby lidar back- 

26 25 

scatter compiled by Russell and Hake and Swissler et al . are of 
interest. From them one can derive the following percentages of time 
that measured peak stratospheric backscatter mixing ratio (x = 0.69 ym) 
occurred in various categories. 
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Percent of Time Measured Values Were: 



Less 

NV 

MV 

Greater 

Time 

Than 

to 

to 

Than 

Span 

NV 

MV 

R18 

R18 

1964-73 26 

5 

70 

13 

13 

1974-80 25 

50 

39 

2 

8 

1964-80 25 * 26 

23 

57 

8 

11 


NV = nominal midlatitude nonvolcanic model (average of 
high- and low-latitude). 

MV = midlatitude moderate volcanic model. 

R18 = Reference 18 model. 

Note that measured values equaled or exceeded the Ref. 18 model value 
only 11 percent of the time in the 17-year period from 1964 through 
1980. For this same period, measured values occurred most frequently 
(57 percent of the time) between our nominal nonvolcanic and moderate 
volcanic values. However, measured values were actually less than our 
nominal nonvolcanic value 23 percent of the time in 1964-80 and 50 per- 
cent of the time in 1974-80. 

When one considers that the above model values are peak values, and 
that the peak is usually confined to a layer several km thick or less, 
it is evident that aerosol particle backscatter contamination in all but 
a few km of the upper troposphere and stratosphere, at wavelength 0.35 
pm, should most of the time be several percent or less, and should very 
rarely attain the values of Ref. 18. In addition, the simultaneous 
availability of the 1.06-pm particle measurements immediately highlights 
the heights and times when particle contamination occurs, and it also 
permits an approximate correction for those effects. 
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